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Uber die Ausbreitung von Zylinder- und Kugelwellen 
in reibungsfreien Gasen und Flüssigkeiten 


Von Karl Bechert 


Zusammenfassung 


Für das Problem der Ausbreitung von Zylinder- und Kugelwellen in 
reibungsfreien Gasen und Flüssigkeiten bei polytroper Zustandsänderung werden 
strenge raumzeitlich veränderliche Lösungen angegeben. — I. Einleitung, 
Ergebnisse. II. Zylinder- und Kugelwellen, bei denen g = f(w). $1. Zurück- 
führung der Aufgabe auf die Integration einer gewöhnlichen Differential- 
gleichung erster Ordnung (16). $2. Allgemeine Eigenschaften der Bewegung. 
Die singulären Integrale von (16). $3. Zur Integration von (16); der Fall n = — 1 
(Flüssigkeiten). $4. Eine spezielle Zylinderwelle jn Flüssigkeiten. III. All- 
gemeinere Zylinder- und Kugelwellen. $5. Ein allgemeinerer Ansatz für 
Zylinder- und Kugelwellen; Zurückführung auf die Integration einer gewöhn- 
lichen Differentialgleichung erster Ordnung (51). § 6. Allgemeine Eigen- 
schaften der Bewegung. Der Fall a =—1. Die singulären Integrale von (51). 
$7. Einiges über die Gleichung (51). IV. $8. Das Problem der allgemeinen 
Lösung für Zylinder und Kugelwellen. — Anhang: Über die Ableitung des 
allgemeinen Integrals aus dem vollständigen Integral (52). 


I. Einleitung, Ergebnisse 


Einleitung. Die vorliegende Arbeit handelt von der Ausbreitung 
von Zylinder- und Kugelwellen in reibungsfreien Gasen und Flüssig- 
keiten; auch die Wärmeleitung soll vernachlässigt werden, ebenso 
die Schwerkraft. Über den Druck machen wir die Voraussetzung, 
daß er nur von der Dichte abhängt, wie es für alle polytropen 
Vorgänge, also auch für adiabatische, zutrifft. Strenge Lösungen 
dieser Aufgabe sind bisher nicht bekannt, soviel ich weiß. Die 
allgemeine Lösung wird in dieser Arbeit auch nicht gegeben werden, 
wohl aber die mathematische Beschreibung einer ganzen Klasse von 
Zylinder- und Kugelwellen, die strengen Lösungen der Aufgabe 
entsprechen. 

Zylinder- und Kugelwellen behandeln wir gemeinsam; r bedeutet 
bei Kugelwellen die Entfernung vom Mittelpunkt der Kugel, bei 
Zylinderwellen die Entfernung von der Zylinderachse; von der 
Strömungsgeschwindigkeit » sei nur die radiale Komponente u von 
Null verschieden. o ist die Dichte, p der Druck; o, u sollen nur 
von r und der Zeit t abhängig sein (Zylinder- bzw. Kugelsymmetrie 
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der Wellen). Partielle Differentialquotienten bezeichnen wir dureh 


ou 
untere Indizes; z. B. u, = er 


Die hydrodynamischen Bewegungsgleichungen und die Konti- 
nuitätsgleichung lauten in beliebigen rechtwinkligen Koordinaten: 


(la) grad bv? — rot + - grad p = 0; 


(1b) + (v grad)o +o = 0, 


Nach unseren Voraussetzungen ist rot ) = 0 und nur die r-Kompo- 
nente der linken Seite von (la) von Null verschieden; sie gibt: 


(2a) u+uu+ - p,= 9; 


aus der Kontinuitätsgleichung (1b) wird: 


ou=0; 


2b) we, + ou, + 
k ist gleich 2 für Zylinderwellen, gleich 3 für Kugelwellen; setzt 
man k= 1, so geben (2a, b) formal die Gleichungen der ebenen 
Wellen, wenn r als kartesische Koordinate gedeutet wird; dieser 
Fall ist in den Ann. d. Phys. [5] 37. 89. 1940 und 38. 1. 1940, 
behandelt. 

Wie in den eben genannten Arbeiten nehmen wir an, daß: 


(2 Cc) p Po = 


daß die Zustandsänderungen also nach dem Polytropengesetz ver- 
laufen; p, ist eine Konstante. Die Adiabaten der idealen Gase 


erhält man aus (2c) fir ,=0 und n=x= ie ; eine brauchbare 


Zustandsgleichung der Flüssigkeiten entspricht p, +0, n =—1; 
n setzen wir mit dem Ansatz (2c) als von Null verschieden voraus. 
Von der Integration der Gleichungen (2a, b) unter der Voraus- 
setzung (2c), und von den Eigenschaften ihrer Lösungen wird in 
dieser Arbeit die Rede sein. 
Es ist zweckmäßig, statt 9, u die Veränderlichen einzuführen 
(n ist ungleich 1 angenommen) 
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(3) E= 


MR 
dani 
(4b) 
zieli 
| die 
dant 
(k - 
is 
Glei 
Die 
ist, 
kons 
schw 
ce ist 
von | 
erhäl 
alle ı 
| 


K. Bechert. Ausbreitung von Zylinder- und Kugelwellen usw. 171 


dann wird aus (2a, b), wenn wir (2c) in (2a) einsetzen: 
kVn-1 


auch aus (3), (4a, b) erhält man die für ebene Wellen gültigen Be- 
ziehungen, wenn k = 1 gesetzt wird? 

Ergebnisse. In Teil II wird vorausgesetzt, daB &=F(n), daß 
die Dichte also nur von der Strömungsgeschwindigkeit abhängt. 
Das Problem der Integration der Differentialgleichungen (4) kann 
dann zurückgeführt werden auf die Integration einer gewöhnlichen 
Differentialgleichung erster Ordnung (§ 1) für die Funktion q von g: 


Ey 
2 r 


= 0: 


l +1 2 
g ist ein Parameter und hängt mit &, 1, r, t zusammen durch die 
Gleichungen (§ 2): 


4 = C exp (JS): E=yı; 


r n kin — 1) 
= 1 
Vn-1 | 2 #41) 


Die Lösung von (4), die in den vorigen Gleichungen dargestellt 
ist, enthält 3 willkürliche Konstanten: C, r und die Integrations- 
konstante der Differentialgleichung für q. 

Vorgegebene Werte von Dichte und Strömungsgeschwindigkeit 
schreiten bei diesen Wellen mit zeitlich unveränderlicher Ge- 
schwindigkeit fort (§ 2); bei adiabatischen Vorgängen in Gasen ist 
diese Geschwindigkeit: 
dé 


dr 
=u-+cyk 


dt 
x—1 

c ist die Schallgeschwindigkeit: c= ag 2 . Je nach dem Vorzeichen 
von d&/dm laufen die Kugeln oder Zylinder vorgegebenen §- und 
n-Wertes schneller oder langsamer als die Strömung u. Für k=1 

erhält man die Formel für ebene Wellen (8 2): 
dr 
dt 
Stoßwellen gibt es beim Zusammenhang £ = F(n) entweder für 

alle t oder gar nicht (§ 2). 


=ute. 


1) Vgl. Ann. d. Phys. [5] 37. S. 94. 1940 GI. (16) und (20). 
12* 


ti. 
n: 
aD 
A 
0, 
re 
1; 
is. 
in 


172 Annalen der Physik. 5. Folge. Band 39. 1941 


Aus der singulären Lösung der Differentialgleichung für g(y) 
ergibt sich das weitere Integral: 
Vn-1 
kn-D) 


9 
- 


7; 


= 
1+ 


es tritt gleichberechtigt neben die oben angegebenen (regulären) 
Lösungen ($ 2). 

Die Differentialgleichung für q läßt sich so umformen ($ 3), 
daß die Integration für n =— 1 in geschlossener Form gelingt; es 
gilt dann: 

(kg? DEI. gr 


const; n=— 1. 


(ka p — 
Ebene Wellen verzerren sich für n = — 1 beim Weiterwandern nicht; 
Zylinder- und Kugelwellen aber wohl (§ 3), natiirlich auch bei n + — 1, 
Für n+—1 ist in §3 eine Form der Differentialgleichung an- 
gegeben, die zur numerischen Integration bequem ist. In § 4 wird 
die Lösung für eine spezielle Zylinderwelle in Flüssigkeiten (n = — 1) 
explizit angeschrieben, und die Möglichkeiten fiir Stoßwellen unter- 
sucht. 
Teil III bringt eine Verallgemeinerung des Ansatzes & = F (n); 
es soll jetzt gelten: 
& = F (re 1), 


Dieser Ansatz liefert bei beliebigem « Lösungen ($ 5). Auch hier 
läßt sich die allgemeine Lösung angeben; sie gibt r«n, r«£, ah 
in Abhängigkeit von einem Parameter » [§ 6, Gl. (52)] und enthält 
4 willkürliche Konstanten; man hat damit also ein vollständiges 
Integral des ursprünglichen Problems (4) gefunden. Dabei ist stets 
«+—1 vorausgesetzt. — Um die Lösung wirklich auswerten zu 
können, muß man eine gewöhnliche Differentialgleichung erster 
Ordnung für eine Funktion q(g) integrieren [$ 5, Gl. (51)]. 

Vorgegebene Werte von r*é und rn schreiten bei diesen Wellen 
mit veränderlicher Geschwindigkeit fort, die von r wie r== und 
von £ wie (( + 7)*+! abhängt (§ 6); r ist eine Konstante. Die Wellen 
verzerren sich im allgemeinen beim Weiterwandern. Stoßwellen 
gibt es entweder für alle ¢ oder gar nicht (§ 6). Neue singuläre 
Integrale liefert der in III betrachtete allgemeinere Ansatz (mit « + 0) 
nicht (§ 6). 

Der Ausnahmefall « =—1 ist in (55), (56), $ 6 behandelt. 
Hier schreiten vorgegebene Werte von „/r und &/r ebenfalls mit 
veränderlicher Geschwindigkeit fort, die proportional r ist und vont 
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exponentiell abhängt. Die Wellenverzerrung zeigt sich wieder. Stoß- 
wellen gibt es wieder entweder für alle t oder gar nicht. — In §7 
steht einiges Wenige über die Differentialgleichung (51). 

Teil IV handelt vom Problem der allgemeinen Lösung für 
Zylinder- und Kugelwellen, und zwar von der Frage: Kann man 
aus dem vollständigen Integral (52) das allgemeine finden? Bei einer 
partiellen Differentialgleichung erster Ordnung für eine unbekannte 
Funktion von zwei Variablen ist die Frage bekanntlich leicht zu 
beantworten; man hat dann zwei willkürliche Konstanten im voll- 
ständigen Integral. Betrachtet man die eine Konstante als Funktion 
der andern und bildet die Einhüllende der so entstehenden ein- 
parametrigen Flächenschar, so gewinnt man durch Elimination des 
Parameters das allgemeine Integral. An einem speziellen Integral 
für ebene Wellen wird gezeigt, daß das Verfahren bei simultanen 
Differentialgleichungen zwar in besonders günstigen Fällen allgemeinere 
Integrale aus speziellen machen kann, daß es aber beim Problem 
der Kugel- und Zylinderwellen nicht brauchbar ist ($ 8). — Im An- 
hang wird eine Methode besprochen, — sie stammt aus dem Lehr- 
buch von Forsyth und ist im Grunde die Methode der Variation 
der Konstanten — die grundsätzlich vom vollständigen Integral zum 
allgemeinen führt. Aber die dabei zu machenden Rechenschritte 


sind in unserem Fall wohl nicht ausführbar, und das System von 
zwei simultanen partiellen Difierentialgleichungen erster Ordnung, 
das am Ende zu lösen wäre, — wenn man es aufstellen kann — 
würde aller Wahrscheinlichkeit nach nicht leichter zu integrieren 
sein als das ursprüngliche System (4). Das Problem der allgemeinen 
Integration für Kugel- und Zylinderwellen bleibt also offen. 


Il. Zylinder- und Kugelwellen, bei denen o nur von u abhängt 


Bei ebenen Wellen bekommt man durch den Ansatz eines Zu- 
sammenhangs zwischen Dichte o und Strömungsgeschwindigkeit u 
sehr allgemeine Lösungen, die eine willkürliche Funktion enthalten. 
Bei Zylinder- und Kugelwellen ist dies nicht der Fall; das hängt 
damit zusammen, daß in (4b) das Glied mit 1/r für k +1 nicht 
verschwindet. Doch ergeben sich auch jetzt noch Lösungen, die 
zwar keine willkürliche Funktion, aber immerhin 3 verfügbare 
Konstanten enthalten. 


$1. Zurückführung der Aufgabe auf die Integration 
einer gewöhnlichen Differentialgleichung erster Ordnung (16) 


Wir setzen also an: 
(5) & = F(n) 


| 
| 
| 
| 
t 
t 
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und erhalten durch Einsetzen in (4a, b): 
nn kVn—-1 dF 
a) Yn-1 +77 
(k— YVn-1 Fy 
2 


re 


Multiplizieren von (da) mit dF/dy und Abziehen von (5b) gibt die 
Gleichung: 


aus der man schließen kann, daß r»,/y nur von 7 abhängt: 

dyn _ dr, 

r+ fit) ist daher Funktion von 7 allein, oder umgekehrt: 

(4) 4 =G(r- fit). 

Bei (6) trennt sich der Fall k = 1 grundsätzlich von den andern 
(k+1); k=1 gibt F =+ 7 + const. 

Da wir bis jetzt aus (da, b) erst eine Gleichung, nämlich (6) ge 
schlossen und ausgenützt haben, so können wir (5a) oder (5b) dazu 
benützen, f(t) genauer zu bestimmen. Einsetzen von (7) in (5a) gibt, 
wenn wir die Abkürzung einführen: 


(8) L=r-fit) 
die Gleichung: 


dG df G u 
9) 


es ist = + 0, sonst wäre „7 raum-zeitlich konstant, folglich muß) 


Funktion von 4 allein sein, weil alle übrigen Glieder in (9) nur voné 

abhängen. Also ist notwendig: 

ft dt 

A und r sind willkürliche Konstanten. Aus (6) wird so: 


dG (a) -1 (k-1)G 


(6a) = g(n), oder 


rn 
n 


1 


(10) =const=— 4; f= STETE 


di dG 


= (0 schließen wir als uninteressant aus; es würde Zeitunabhängig- 


keit von n und & bedeuten. 
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und aus (9): Yeni 
kVu —1 dF 
Po 
(11) n 2 dG 
7 laBt sich eliminieren; nach (11a) ist: 
dG \* di, 
= 
und nach (11b): 
km-—1) ad dF 
(# ie) 16 
2 dG 
Gleichsetzen gibt eine Beziehung für F in Abhängigkeit von G: 
d’F adF \? n+1 
| (k — 1)GF dG? +6 (de) 
dF \* dF 2G 


n—| 


(12) 


die Gleichung ist homogen in F und G, daher ist der Ansatz zweck- 
mäßig: 
(13) Feaq-G; G=e; 


so kommt: 


(14) 
do 


Wir führen gy als unabhängige Variable ein und 
(15) 
als abhängige, dann erhalten wir endgültig: 
(16) 9) +4 (k— + =0, 
eine gewöhnliche Differentialgleichung erster Ordnung fiir q. 
$2. Allgemeine Eigenschaften der Bewegung. 


Die singulären Integrale von (16) 


‚Wenn die Gl. (16) gelöst ist, also q in Abhängigkeit von ¢ 
bekannt ist, sind auch 7, 2, € in Abhängigkeit von y angebbar. Es 
folgt nämlich aus (15), (13) und (7): 


(17a) 4 =G =C exp (f 
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C ist eine beliebige Konstante. Aus (13) und (5) kommt: 


(17b) 
und aus (11b), (13), (15), (8), (10): 

7 r 1) 
Wir haben eine Darstellung der Dichte wa der Strémungsgeschwin- 
digkeit (die aus &, 1 gemäß (3) zu finden sind), und der Raum-Zeit- 
Abhängigkeit Fr durch den Parameter g. Wenn also (16) gelöst 


ist, kann man &, n, punktweise berechnen, indem man ¢ 


t+r 
punktweise wählt. 

Unter der von uns gemachten Voraussetzung, daß die Dichte 
nur von der Strömungsgeschwindigkeit abhängt, geben die Gl. (17) 
die allgemeine Lösung des Problems; sie enthält 3 willkürliche 
Konstanten, nämlich C, r, und die Integrationskonstante B [vgl. 
Gl. (30a, bj], die bei der Lösung von (16) auftritt; die Konstante A 
der Funktion f(t) von G1.(10) ist herausgefallen. Setzt man die Werte 
der Konstanten fest, so ist damit &, 4 für t=0 als Funktion 
von r gegeben; (17c) bestimmt dann das Verhalten für spätere t. 

Aus der Annahme (5): $=F(n) folgt, daB § = rei zugleich 


n = const bedeutet, und aus (7), (10), daß dann auch konstant 


ist. Feste (d. h. bestimmt, aber beliebig, gewählte) Werte “ie Dichte 
und der Strömungsgeschwindigkeit bewegen sich also bei diesen Wellen 
mit konstanter Geschwindigkeit dr/dt weiter; dr/dt bestimmt sich 
aus (17c) zu [vgl. (3)]: 


(18) 


dabei ist zu beachten, daß auch @ konstant ist und damit q, wenn 7 
oder & konstant sein soll. 

Die Gl. (18) kann man übrigens auch direkt aus (4a) finden, ohne 
Kenntnis des Integrals (17). Dazu hat man nur zu beachten, dab 


aus 7 = const folgt: 1, + 7, = 0 und daß 


ist. (18) läßt sich anschaulich deuten. Der Ausdruck u - ¢4 
k Vn—-1 


ist nämlich gleich: + 


S , also gilt wegen (3): 


1 d 
(18a) utag ? yh 
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bei adiabatischen Vorgängen in Gasen bedeutet der Faktor vor Yk = 


die Schallgeschwindigkeit c. Für adiabatische Vorgänge in Gasen 
ist daher die Geschwindigkeit des Fortschreitens fester Werte der 
Dichte und der Strömungsgeschwindigkeit gegeben durch: 


(18b) 


Je nach dem Vorzeichen von dé/dy laufen die Kugeln oder Zylinder 
konstanter £- und 7-Werte schneller oder langsamer als die momentane 
Strömung u. Der Spezialfall ebener Wellen ergibt sich aus (18a) 


oder (18b), wenn k = 1 gesetzt wird; dann wird nämlich In =+1, 


— vgl. die Bemerkung nach Gl.(7) —, und es kommt: an =ute 


Ohne ausführliche Diskussion kann man an (17c) nicht er- 
kennen, ob & = mv bei festem ¢ eine monotone Funktion von r wird 
oder nicht. Wenn £ nicht monoton ist, gibt es bei festem ¢ und & 
zwei oder mehr r-Werte, welche die Gl. (17c) erfüllen; die Funktion 


= (14+ 1) $4) 


auf der rechten Seite von (17c) ist ton notwendig verschieden für 
diese verschiedenen r-Werte. Daher ist auch die Geschwindig- 
keit dr/dt der zugehörigen &-Werte verschieden; jedes einzelne dr/dt 
ist für sich konstant, nach unseren früheren Überlegungen. Das 
bedeutet: Wenn es nicht-monotone Funktionen £(r) zu festem t 
gibt, dann verändert sich ihre Breite im Laufe der Bewegung pro- 
portional der Zeit. Unter Breite ist dabei der Abstand zwischen 
zwei Punkten gleichen &-Wertes bei festem ¢ verstanden. Die Pro- 
portionalität mit der Zeit ist die unmittelbare Folge der konstanten 
und verschiedenen dr/dt dieser Punkte. 

Im allgemeinen wird die „-Kurve und daher auch die &-Kurve 
ihre Form bei der Bewegung ändern, denn dr/dt ist von » abhängig. 

Eine Stoßwelle bekommen wir, wenn 7,, 9,, &,, & unendlich 
werden. Wir bezeichnen die rechte Seite von (17c) mit ®(g) und 
berechnen: 


= n,= = 
Or dy Ir ®yn-1i du 


nach (18). Andererseits ist: 


dig®@ Me (=) Nr 


dt 


kn-l) ad dy 
dp 


re 1) + a) — 


a dé 
= 
| 
d (dr k(n —1) 
kin-1 
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nach (17a). Einsetzen von (16) gibt schlieBlich: 

d Ir l e kin — 1) 
Daher: 
(18d) = 


Ebenso: 


nk—-D 
r(kq?—1) ’ 


ot 
Das gibt: 


n?(k — 1) (1 + 


18e) 7, =— 


Aus (5), (7), (13) und (15) folgt: 


(18f) 4 4,- 
Fiir: 


1 
(18g) 


werden im allgemeinen alle vier Differentialquotienten unendlich; 
(18g) ist also eine Stoßwellenbedingung. Bemerkenswert ist, daß 


d (dr be 
du (47) me Vi 


im allgemeinen Null wird; vgl. (18c). 

Die Stoßwellen (18g) gibt es entweder für alle t oder gar nicht. 
Denn diese Gleichung bedeutet die Auswahl eines bestimmten 
q-Wertes, daher eine Festlegung von g, wenn die Integrations- 
konstante von (16) bestimmt gewählt ist; mit gm ist aber 7 festgelegt, 
wenn noch die in (17a) vorkommende Konstante C bestimmt gewählt 
ist. Setzt man schließlich noch z fest, so ist bei vorgegebenem t 
der r-Wert aus (17c) berechenbar, für den die Stoßwelle auftritt; 
da (17c) für alle ¢ gilt, so gibt es immer einen solchen r-Wert, wenn 
die rechte Seite von (17c) reell ist. Ist sie nicht reell, so gibt es 
keine Stoßwellen. 

Mit Hilfe der Ausdrücke (18d, e, f) ist übrigens leicht zu be- 
stätigen, daß die Lösung (17) den Differentialgleichungen (4) genügt. 

Die singulären Lösungen von (16). Sie sind durch das Null- 
werden des Faktors von dgq/dy gegeben. gy = 0 ist uninteressant, 
weil es zu 7 = const, &= 0 führt. Es bleibt die Möglichkeit: 


(19) Ts 
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also nach (15) = = 0; » ist konstant. Den Wert der Konstanten 


erhalten wir durch Einsetzen von (19) in (16); es ergibt sich eine 
biquadratische Gleichung für 4, deren Lösungen sind: 


(20) y=+l; + 


Das zweite Paar Wurzeln ist für n > 1 unbrauchbar, weil g= : 
nach (3) reell sein muß. Fürn =— 1 und k = 2 wird dieses Wurzel- 


paar in $ 4 besprochen. 
Es bleibt g =+1 zu untersuchen. Man bekommt aus (13) 


und (11b): 
(20a) F=H+G; (24+ kin—1)); 
auch diese Lösung ist unbrauchbar, wenn 


(21) n=1—7; 
bei Zylinderwellen würde das n = 0 bedeuten, was wir ohnehin aus- 
geschlossen haben, bei Kugelwellen dagegen n = — — Wenn 


nt1— : findet man durch Einsetzen von (8) und (10): 


n—1 r 

k(n — 1) t+ı = 
i+ 


(22a) y= 


das bedeutet gemäß (3): 
2 


2 mas o=| (n r \n-ı 
2+km-]) t+?" 


a(2+km—) t+r 
Diese Lösung wird man wohl nur zur Schematisierung ganz spezieller 
Strömungsverhältnisse brauchen können. Wenn z.B. für t=0 in 
einem bestimmten Gebiet annähernd gilt: 

(n k ar)" 
und außerhalb wu = 0, eg = 0, so wird man den weiteren Verlauf von 
Strömung und Dichte in diesem Gebiet — dessen Bereich sich 
übrigens zeitlich verändert — mit den Gl. (22) annähernd be- 
schreiben können. 


(22 b) ar, o= (+ 


$3. Zur Integration von (16). Der Fall » = — 1 (Flüssigkeiten) 
Eine spezielle Lösung von (16) findet man durch den Ansatz: 
q=const-g* in der Form: 
2 


9 
=) kin-l)q 
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Zur Vereinfachung der Koeffizienten in (16) wird es zweckmäßig 
sein, q, an Stelle von y als unabhängige Veränderliche einzuführen, 
Außerdem setzen wir: 

(24) W 

was durch die teilweise Homogenität der Gleichung nahegelegt wird; 


die ebenfalls naheliegende Substitution q = = würde eine ver- 
1 


wickeltere Gleichung geben. So 


(25) sen + = (w — 1)(q,? w?— 1). 


Wir vertauschen die nh der Veränderlichen und betrachten q, als 
abhängig. Außerdem setzen wir: 

9 1 
(26) 


und suchen durch geeignete Wahl von z die Gleichung für v einfach 
zu machen. Zunächst erhalten wir: 


Die Forderung: 
(27a) 

führt auf: 

(27b) 

und: 


(k -— D(?2w + D)v 
(270) kww-— 1) 


~ 


Allgemein habe ich diese Gleichung nicht integrieren können. 
Man kann sie noch in mannigfache Formen bringen, von denen einige 
recht merkwiirdige Transformationseigenschaften haben. Numerische 
Integration wird man am besten an (27c) versuchen. Die Isoklinen 
der Gleichung sind fir v als Funktion von w leicht zu zeichnen; 
so ist eine graphische Integration möglich, von der in der folgenden 
Arbeit berichtet werden soll. 

Für einen speziellen Fall, wenn nämlich n = — 1, läßt sich das 
allgemeine Integral von (27c) auch explizit angeben. Wir formen 
dazu (27c) allgemein um durch den Ansatz: 

(28) 


’ 


und bekommen, wenn wir (27a) beniitzen, die Gleichung: 


K. B 


Der ; 
gegeb 
(29a) 
Went 
es wi 
(29b) 
Das : 
(29) 
das 9 
(30) 


Das 


= ( 
erwäl 
stand 
Lösuı 


(30a) 


B ei 
welch 
Die | 


zu tu 
L 


solch 
welle: 
abhäı 
hatte 
verze 
wert, 
zerre! 
gesch 


(31a) 
daraı 
(31b) 


w? dz dv 
(w — 1) dz 4, k-1 
‘s46- & 
2(k—1) 
z=(w—1)* , 
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Der Zusammenhang von U,w mit q, p ist gemäß (23), (24), (26), (28) 


egeben durch: 
k(n — 1) 


1 
(29a) 5 1 


Wenn n =—1, lassen sich in (29) die Veränderlichen trennen und 
es wird: 
9) kw(w — 1)(U—1) 4 200 n=-1. 


Das allgemeine Integral von (29b) ist: 


U(U — 


(w — =const; »=— 1; 


(29 c) 


das gibt nach Umrechnung auf q, y: 
(kg? — g2% 
— 
Das Integral (30) ist auch fiir k= 1 noch brauchbar; dann wird 
g = const = - Daß die Konstante den Wert + 1 hat, ist bei (7) 


erwähnt worden, Für Zylinderwellen in Flüssigkeiten, die der Zu- 
standsgleichung (2c) mit n = — 1 genügen, ist demnach die allgemeine 
Lösung von (16): 
(30a) =B; n=-1l; k=2; 
B eine willkürliche Konstante. Für Kugelwellen in Flüssigkeiten, 
welche der Zustandsgleichung (2c) mit n=— 1 genügen, gilt: 

(3q? — 1) 
(30b) n=-1; k=3. 
Die Konstanten B in (30a, b) haben natürlich nichts miteinander 
zu tun. 

Zusammen mit (17a, b, c) haben wir in (30a) bzw. (30b) für 
solche Flüssigkeiten die Beschreibung aller Zylinder- bzw. Kugel- 
wellen, bei denen die Dichte nur von der Strömungsgeschwindigkeit 
abhängt. Ein Beispiel wird in § 4 gegeben. 

Die Betrachtungen, die wir im Anschluß an Gl. (17) angestellt 
hatten, gelten auch fir n=— 1. Es gibt also auch hier Wellen- 
verzerrung bei Zylinder- und Kugelwellen; das ist darum beachtens- 
wert, weil sich im ebenen Fall für n =—1 die Wellen nicht ver- 
zerren. Die Gleichung (18a) kann nämlich allgemein in der Form 
geschrieben werden: 


(31a) 


daraus folgt: 
(31b) 46-9 @. 


dy?’ 


(30) = const; »=—1. 


| 
dr kYyn-ı 


182 Annalen der Physik. 5. Folge. Band 39. 1941 


bei ebenen Wellen mit & = F(n) gilt aber: § =+ + const. Also 
wird: 


vgl. auch Gl. (29a, b) in der Arbeit I über ebene Wellen'). Die 
Geschwindigkeit des Fortschreitens vorgegebener Werte von wird 
folglich bei ebenen Wellen von 7 unabhängig, wenn n=— 1; d.h. 
alle Werte von n haben dieselbe Geschwindigkeit, die Welle läuft 
unverzerrt fort. 


$4. Eine spezielle Zylinderwelle in Flüssigkeiten 


Wir spezialisieren die Formeln des vorigen Paragraphen so, daß wir 
eine Zylinderwelle in Flüssigkeiten bekommen (n = — 1, k= 2). Dam 
gilt Gleichung (30a), aus der man q berechnet zu: 

1 < 
die Konstante B sei der Einfachheit halber gleich 1 gesetzt, dann 


laBt sich das Integral J ar leicht elementar auswerten. Es wird: 


(33) dg + V)dp 
v2 
Die Vorzeichen in (33) sind einander zugeordnet, man hat entweder 
stets das obere oder stets das untere zu nehmen. Einsetzen von (32), (33) 
in (17a, b, c) gibt schließlich, wenn wir (3) benützen: 


(34) 


D ist eine willkürliche Konstante, welche für die beiden Vorzeichen 
von einander unabhängige Werte hat; außerdem ist zu beachten, 
daB eine Lösung (34) mit bestimmtem D nur entweder für |g! >! 
oder für p|< 1 gilt. Die Funktionen (34) im anderen Gebiet (|g <1 
bzw. g > 1) gehören nicht zum gleichen D wie die vorigen. Das 


hängt damit zusammen, daß in (34) gemäß der Herleitung bei st 


zunächst Klammern an Stelle der Absolutstriche stehen, daß daher 


1) Ann. d. Phys. [5] 37. 99. 1940. 


K 
fü 
m: 
ge 
bil 
(3: 
M: 
(36 
u=D v2, j;0= @ v2, &=1: 
| r jt y2 | 
= = D V2 T 
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eine Lösung mit bestimmtem D nur entweder fir gy >1 oder 
fir gy <1 reell ist. Das Vorzeichen des Ausdrucks von g ist 
mathematisch willkürlich und den übrigen +-Zeichen nicht zu- 


geordnet. 
Die Ausdrücke für u,, u,, o,, o, kann man mit Hilfe von (18d, e, f) 


bilden; dabei gilt wegen B = 1 gemäß (30a): 
(35) 


Man erhält so: 


"yz V2gt + 
; 
ae (g?+@V2+1). 
0,=- — : 
2y2r 


\ 


und für die Abhängigkeit der Geschwindigkeit dr/dt von u nach (18c): 


1 8 
35) 
V2) 
Für das Auftreten von StoBwellen gibt es mehrere Möglichkeiten; 


(6b 


die Forderung (18g) liefert wegen (35): entweder q = g = + = ) 
oder gy =+co; außer (18g) gibt es noch die Möglichkeit: g = + 1. 
Die erste Méglichkeit q = g = + = bedeutet, daß u,, u,, 0,. 0, > © 


gehen, und gleichzeitig u —> 00,0 0, gegen einen end- 


+1 
lichen Wert. Dies ist zugleich der Fall des zweiten Wurzelpaares, 
das in Gl. (20) angegeben ist; diese StoBwelle verhält sich also wie 
eine singuläre Lösung. Die zweite Möglichket g=+w gibt 
wieder U,, 0,, 0, > 00; jetzt geht aber u o und r 


gegen endliche Werte. 


1) Die Vorzeichen hier sind den Vorzeichen in (32) zugeordnet, bei den 
folgenden Möglichkeiten aber nicht. 


ft 
[ + 1 2 2 +V2 4 | 
D V2. (tV2) TP 7 +1 
(ox) 
~ y2 
| | | 
(36a) 
d: | 
er 
3) 
| 
f 
l 
1 
as 
1 
1 
er 
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Die dritte Möglichkeit g = + 1 gibt u, u, > 00, 0, —>0 oder 
umgekehrt, je nach der Verssichenwahl; o, bleibt endlich. Die 
Bedeutung dieses Falles übersieht man am besten, wenn man nach (34) 
allgemein bildet: 


ae (p+ 


t+ı 
: 
(36c) 
für g = + 1 erhält man aus (34) u >» 00 oder 0, 9 ->0 oder o, 
r 


t+. 


Ausdrücke u: Try und o- 437 bei g=+1 endlich bleiben, 


Man bekommt: 
(36d) u=— 


> 20 oder 0, je nach den Vorzeichen. (36c) zeigt aber, daß die 


r as(t + 

was genau mit (22) übereinstimmt; es ist k=2,n=— 1. Inder 
Tat ist ja g =+1 das singuläre Integral, das dem ersten Wurzel- 
paar in (20) entspricht; denn aus (32) folgt, daß hier g= wird. 
Daher verhält sich auch diese (ausgeartete) Stoßwelle wie ein singuläres 
Integral, und zwar so, wie (36d) sich bei t=— r oder bit=+w 
verhält. Eine ausführlichere Untersuchung der Wellen soll in einer 
folgenden Arbeit gegeben werden. 


III. Allgemeinere Zylinder- und Kugelwellen 


85. Ein allgemeinerer Ansatz für Zylinder- und vce 
Zurückführung auf die Integration 
einer gewöhnlichen Differentialgleichung erster Selen (51). 


Als Verallgemeinerung von (5) setzen wir an: 


(37) E= F(g, (r) + m). 

Eine vollständige Diskussion der Möglichkeiten dieses Ansatzes habe 
ich nicht gemacht. Es scheint, daß er ziemlich zwangsläufig auf gy, =r", 
, = 7** führt, wobei « willkürlich bleibt. Im folgenden setzen 
wir daher: 


(38) « beliebig. 
Es läßt sich in ähnlicher Weise wie für «= 0 zeigen, daß dieser 
Ansatz mit den Gleichungen (4a, b) verträglich ist, also wirklich 
Lösungen gibt. 

Wir führen die Abkürzung: 
(39) s-r“.n 


| K. Be 

ein; 
in (48 

(40 a) 

(40b) 

| 
und : 

(41a) 

Also 
wegel 

(41b) 

Wir 

(41e) 
wir I 

(42) 

Wie 

(43) 
dahe: 
(43a) 
Wen 
(43b) 
An 
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ein; es ist gern; = + n,) . Einsetzen von (38) 
in (4a, b) gibt: 


=.) kyn L -a,F 
408) Vn-1 


(k-1VYn-1 


(40b) 
+ Fn, + 9 


Wir multiplizieren die zweite Gleichung mit r“, die erste mit — r“-F, 
und addieren; wenn wir noch (39) beachten, erhalten wir: 


(la) 
Vn -1 kyYn-1 
Also ist notwendigerweise rn,/n Funktion von s allein; d. h.: 
Oilgs—algr) ölgs 


ölgn 


wegen (39). Daraus folgt ähnlich wie in (7): 
(41b) s=G(r-fi). 
Wir benützen wieder die Abkürzung 
(41e) A=r-fi(t) 
und bekommen durch Einsetzen von 7 =r-*-G(d) in (40a), wenn 
wir mit r?«+! multiplizieren und beachten, daß r = - 
4 1 „8 (2. 
a4  Wa-1 


kYn-1 [. dF dG 
+= 4% _ = 


(42) 


Wie in (10) schließen wir, daß f-*-?- ce konstant sein muß: 
(43) 

daher: 

Wenn « =— 1, tritt an die Stelle von (43a): 


(43b) f=e-4+9; u=-—1. 
Annalen der Physik. 5. Folge. 39. 
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Im folgenden setzen wir @#— 1 voraus; von «=—1 wird 
im nächsten Paragraphen die Rede sein. 
Aus (42) wird so: 


+2 da G 
dh Va-1 


| (az dG 


| Ais 
(44a) 


dG =0. 


Aus (41a) entsteht wegen —a nach Multiplikation 


mit — =! die Gleichung: 


(k—1)(n—1)—- 2a ‚ar k(n — 1) dF \* 
2 P+ - 
(n — 1) dG dF \? 
kn-1l) dF 
(44a, b) sind zwei gekoppelte gewöhnliche Differentialgleichungen 
erster Ordnung für F,G als Funktionen von A, denn F hängt 
über G von A ab. « ist willkürlich geblieben, wie wir behauptet 
hatten; der Ansatz (38) gibt für beliebiges « Lösungen. Wenn « = 0, 
folgen aus (44a,b) die Gl. (11b,a), wie es sein muß. Die Gl. (44a,b) 
erlauben die Transformation: 


Dia, G—> Deti.G, F—> Deti.F, 
mit der beliebigen Konstanten D. Daher ist es zweckmäßig, als 
neue Veränderliche: 
(45) u 
di 


— 


dF a F du 
| +- ig) 


einzuführen. Multiplizieren von (44a) mit - Yn — 1 gibt: 


(46a) 


G di 


iF dG and bekommen: 


(44b) multiplizieren wir mit 


2 tes a *G 


G du dF \? Pr 
| “at Du dG (5G) ) =o 


| " — 1)(in — 1) —a(n +1) G dF kin — 1) F aF 
(46 b) 


K. 
Da 
chi 
(47 
di 
(47 
| 
| 
(48) 
auf 
(49) 
und 
(50: 
Das 
(501 
Wir 
ag 
du 
| Ay \ (a + 26) q al 
(b0¢ 
| und 
(50d 


K. Bechert. Ausbreitung von Zylinder- und Kugelwellen usw. 187 


Damit ist @ als unabhängige Veränderliche eingeführt. Die Glei- 
chungen sind jetzt homogen in F,G, u; also wird man weiter setzen: 


(47) F=q-G; p=y-G; G=e. 
Das gibt: 


(47a 
| 


y do 


e+ 


(k —1)(r— 1) —a(n4+1) “ dq kin — 1) 
| 


47b) 


y do 


y kann eliminiert werden. Wir führen wie in (15) ein: 
(48) g= +9, 


außerdem die Abkürzungen: 


und erhalten aus (47b): 


=): 


(50a) PR dig -1+ N 


— 


Das setzen wir in (47a) ein und bekommen: 


Z 


(50b) 
| Z= (N — (1 

Wir fassen jetzt g als unabhängige Veränderliche auf und bilden: 
“r_ == i (q—g) mit Hilfe von (50b). Einsetzen dieses Aus- 


drucks in (60s) gibt eine Differentialgleichung erster Ordnung fiir 
q als Funktion von g. Der Reihe nach kommt: 


(0) 
und als Differentialgleichung: 
dZ dN + 
(60d) @— +Z(N + —kq’) =0. 
i3” 
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Weiter: 
dZ aN, di ay 
Fr +4) +7, (1 — kg’) (1 + 37g) 


gy 


20744; 


dq 
NdZ _ ZdN . 


der Faktor von m = ist: 
dp dp 


a? 
(50e) (3° p— (l+y7 
Die Differentialgleichung wird so: 


dı » a? 2 


(1) (Na+ + 879%) 


für «= 0 geht diese Gleichung in (16) über, wie es sein muß. 


$6. Allgemeine Eigenschaften der Bewegung. Der Fall « = - 1. 
Die singuliren Integrale von (51) 
Wenn (51) gelöst ist, findet man &, 7, A in ähnlicher Weise 
wie in $2. Aus (39), (41b), (47), (48) folgt: 
(52a) n = 1r-*-C-exp (/ 
aus (38), (47): 
(52b) 9-4; 
und aus (41c), (43a), (45), (47), (50b), (50c): 
| ttt Vn-1 N 
_ (@+)) + yo’) 
| Vn-1 “| kN 
Der Ausdruck N ist aus (50a) noch einzusetzen; er ist als Funktion 
von p bekannt, wenn (51) integriert ist. Für «= 0 gehen die 
Gl. (52) in (17) über. 
Wir haben wieder, wie in $2, eine Darstellung der Dichte, der 
Strömungsgeschwindigkeit und des raum-zeitlichen Verhaltens durch 


den Parameter 9. Diesmal enthält die Lösung aber vier willkürliche 
Konstanten: «, C, r und die Integrationskonstante von (51), die 


dp 


(52c) 
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Li 
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in q steckt; A ist wieder herausgefallen. — (52) gibt die allgemeine 
Lösung der Aufgabe (4) unter der Voraussetzung (38). 

Konstantes r*-& bedeutet hier [vgl. (38)] zugleich konstantes 
r«.7 und daher nach (52b) konstantes gw; dann bleibt nach (52c) 


a+l 
konstant. Feste Werte von r°-=B, und von w-y=b, 
T 
schreiten also bei den Wellen der Form (38) mit der Geschwindig- 
keit fort: 


| dr Zn Z B, 


a 


dt N 
1 
a+l 


(53) 
BZ 


sie ist nicht konstant, hängt vielmehr von r wie r~* ab, oder, weil 


konstant ist, von ¢ wie + 


Da Zn/N im allgemeinen nicht von »; unabhängig ist, so hängt 
dr/dt von n ab; die Wellen verzerren sich also beim Weiterwandern. 
Denkt man sich r*& bei festem t als Funktion von r++! auf- 
getragen, so kann r@& monoton sein oder nicht (vgl. die entsprechen- 
den Betrachtungen in $ 2). Ist r*& nicht monoton, so ist nach 
(52c)Z/N g mehrdeutig. Gemäß (53) ist die Wanderungsgeschwindig- 


a+i 


a+1 
keit eines festen Wertes von r*é zeitlich unveränderlich. 


Die Breite der Kurven r«£ (in Funktion von r*+?!) ändert sich 
also proportional der Zeit, wenn r*& eine nicht-monotone Funktion 
von ist. 

Die Berechnung der Differentialquotienten 7,, 7, &.&, wird im 
Anhang unter der Voraussetzung durchgeführt, daß auch die vier 
Größen, die wir hier als Konstante erhalten haben [«, r, C und die 
Integralkonstante B von (51)], von r,t abhängig sind. Die Formeln, 
die wir jetzt brauchen, gewinnen wir daraus, indem wir die Ab- 
leitungen von «, r, C, B gleich Null setzen. So kommt mit der 
Abkürzung: 

Vn-1 N 
(688) oi 
aus den Gl. (f’, f”), ig‘, g”) des Anhangs: 


“(1-0 


&=40,; §& =90,+ 


| 
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Wegen (50d) wird: 
dig: _ _ — kq”) 
l (q f) “dg N ’ 
und schließlich: 
2n N is. 


(n r Vn- 1’ or (kg? -1)’ 


2nNq 2Z 474 
Va - 1° ork?—-) 


(53c) 


Für «= 0 erhält man aus diesen Gleichungen wieder (18d,e, fi. — 
Daß durch (53c) die Differentialgleichungen (4) identisch erfüllt 
werden, ist ohne große Mühe zu bestätigen. 

Auch hier ist offenbar 


(53d) q=+ Vi 


im allgemeinen eine StoBwellenbedingung, und die zugehörigen 
Wellen gibt es entweder für alle ¢ oder gar nicht; das kann ebenso 
wie bei Gl. (18g) erkannt werden. 

Der Fall «=—1. Bei allen Formeln dieses Paragraphen ist bis 
jetzt — 1 vorausgesetzt gewesen; «= — 1 war von (43) an aus- 
geschlossen. Die Gl. (44a, b) sind aber auch für «= — 1 noch richtig. 
Die Transformationseigenschaft von (44) ist aber jetzt: 


6—6, F—+F; 
sie legt nahe, zu setzen: 
= u; 


dann folgt aus (44a) nach Multiplikation mit u Yn — 1: 


(54a) —Ayn—1 +G(14+G te) + (3 + F =0: 


G dai. 
Fi dG’ 


i dF F\? 1 - dF 
4 ed (1 +7(@))|@ + (n 1) (1 — k (46) 


du/dG kann leicht eliminiert werden und es ergibt sich eine Difie- 
rentialgleichung erster Ordnung für F in Abhängigkeit von G: 


| G ig (1 +7(@))+ 


und aus (44b) nach Multiplikation mit 


(55) 
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Wenn (55) integriert ist, erhält man die endgültige Lösung in 
ähnlicher Weise wie vorher. Es wird jetzt: 
(56a) n=r-G; 
(56b) £E=r.F; 
aus (54a) oder (54b) ist du/dG berechenbar, also auch u und daher 2. 
Das gibt gemäß (41c) und (43b): 
G 
r-exp(— At)=C- exp dG}. 


(e-4: liBt sich mit C zu einer Konstante vereinigen). Die Gl. (56) 
geben eine Darstellung der physikalischen GréBen durch den Para- 
meter G. Die Lösung enthält nur 3 Konstanten: A,C und die 
Integrationskonstante B von (55). (56) ist die allgemeine Lösung 


unter der Voraussetzung: §=r-F (2 


Feste Werte von n/r und von &/r schreiten hier mit der Ge- 
schwindigkeit fort: 


auch hier ist die Geschwindigkeit nicht konstant; sie hängt von r linear, 
daher von ¢ exponentiell ab. Auch hier hat man im allgemeinen 
Verzerrung der Wellenform beim Weiterwandern. 

Gibt es zu festem &/r (d. h. F) und zu festem t mehrere r-Werte, 
so ist &/r eine nicht-monotone Funktion von r, daher auch von lg r. 
Nach (57) ist as A 
ändert sich die Breite der Kurven &/r (als Funktion von lgr) pro- 
portional der Zeit, wenn &/r eine nicht-monotone Funktion von lg r ist. 

Die Differentialquotienten 7,, 7,, £,, &, werden fiir @ =— 1: 
dF 
dG 


AVn-1-G-y7F 


für festes &/r zeitlich unveränderlich; folglich 


AYn-1G-yFG,,+yF° 
dG 


dF dF 
AYn IF r A(G 


R 
| Ayn-1-G—-7F AYn-1-G-7F 


Das Nullwerden des gemeinsamen Nenners liefert eine Stoßwellen- 
bedingung; auch hier gibt es die zugehörigen Stoßwellen entweder 
für alle t oder gar nicht. — Daß die Ausdrücke (57a) die Gl. (4) 
erfüllen, erkennt man durch Einsetzen, wobei (55) zu beachten ist. 


rA(G? +yF3ı 


N, = 1,=— 


t 


(57a) 4 


G dF 
16); 
(7) A edt AC exp 
| 
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Die singulären Lösungen von (55) bekommt man durch Null 
setzen des Faktors von dF/dG. G = 0 ist unbrauchbar; Nullsetzen 
der eckigen Klammer liefert F als Funktion von G; damit soll das von 
dF/dG freie Glied in (55) ebenfalls Null werden. Das gibt G = const, 
F = const; also müßte sein [vgl. (56a, b)]: 7=r - const; =r - const, 
Einsetzen in (4) zeigt aber, daß in Wirklichkeit keine solchen Lösungen 
existieren. 

Die singulären Lösungen von (51). Für die singulären Lösungen 
von (51) muß der Faktor von dq/dg in dieser Gleichung verschwinden, 


Das gibt zunächst die Möglichkeiten: l.q=@; 2.gy=0;3. + an = 


1+y7g’. Einsetzen in (51) zeigt, daß 1. bedeutet: entweder q=q =0 


— das heißt z = 0, was wegen (38), (39), (41b) eine für uns uninter- 


essante Lösung ist — oderg=y=+ y | __ das führt gemäß (47a) 


_ 


und (48) auf y = 0, d. h. nach (47), (45) A = 0, was auch uninteressant 
ist — oder q = y= +1; 2. gibt wieder a = 0; 3. bedeutet p = const, ' 


daher nach (51) auch q = const und weiter — vgl. (48) — tw const; 


do 

die Konstante muß Null sein, weil sonst g nicht konstant wäre. 
Daraus folgt aber weiter, dab q = ist und das ist der Fall 1. 
Es bleibt also im ganzen nur die Möglichkeit q = g =+ 1 übrig 
und diese führt durch Einsetzen in (47a, b) wegen « +—1 auf 

dw 

=0; ya 

do Atfe+1)VYn-1 
und schließlich genau auf die Gl. (22), also auf dieselbe singuläre 
Lösung wie im Fall «@=0; «+0 gibt demnach keine neuen 
singulären Lösungen. 


87. Einiges über die Gl. (51) 


Die Gl. (51) hat, wie schon (16), ein leicht angebbares partikuläres 
Integral: Z= 0, das für « = 0 in (23) übergeht. Nach (50d) kann 


nämlich mit Z = 0 zugleich = = 0 erfüllt werden; dann folgt aus 


(50d) weiter, das alle Differentialquotienten von Z nach @ Null sind. 
Z = 0 ist daher für alle g mit (50d) verträglich, also ein Integral 
von (50d). Es lautet nach (50c): 
y(k— 1) 
ER 
(58) -— 2 ) 
aty 


k 
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Physikalisch brauchbar ist dieses Integral leider nicht, weil es 


wegen (52c) ; 0 bedeutet. 
Wenn n =—1, läßt sich (51) vermutlich ähnlich wie (16) in 
geschlossener Form integrieren. Man wird zweckmäßig den Ansatz 


machen [vgl. (30)]: 
9 kg? —1)-(4)"- [a - g?— He, g) = B: 
59) (kg’ — 1) ¥) 


ich habe die Funktion H aber nicht berechnen können. Im allgemeinen 
Fall (n #— 1) wird man (51) numerisch integrieren müssen. 


IV 
§ 8. Das Problem der allgemeinen Lösung 
fiir Zylinder- und Kugelwellen 

Das Integral (52) enthält 4 verfügbare Konstanten: «, C, r und 
die Integrationskonstante B von (51), ist daher ein vollständiges 
Integral des Gleichungssystems (4). Es liegt nahe zu versuchen, aus 
ihm das allgemeine Integral in ähnlicher Weise aufzubauen, wie dies 
bei den partiellen Differentialgleichungen erster Ordnung mit nur 
einer unbekannten Funktion möglich ist. Bei den ebenen Wellen ist 
dieses Programm in der Tat wenigstens soweit leicht durchführbar, 
daß man Lösungen mit einer willkürlichen Funktion bekommt. Man 
benutzt dazu die singuläre Lösung (22a); zunächst bestätigt man 
leicht, daß (22a) auch noch für k = 1 Lösung von (4) ist; sie enthält 
nur eine verfügbare Konstante. Eine weitere Konstante gewinnt man 
aus der Bemerkung, daß k = 1 das Glied mit 1/r in (4b) zum Ver- 
schwinden bringt; dann erlauben die Gl. (4) die Substitution (statt r 
schreiben wir ©): 2 —> 2 + const, wie ja physikalisch unmittelbar 
einleuchtet; es ist kein z-Wert mehr vor dem andern ausgezeichnet. 
Daher ist auch: 
2Vn—-1l 
Lösung von (4), wenn & eine beliebige Konstante bedeutet. Eine 
dritte Konstante ergibt sich daraus, daß bei ebenen Wellen der all- 
gemeine Ansatz: &= + 7+ const auf Lösungen führt‘); addiert man 
zu den rechten Seiten von (60a) zwei zunächst beliebige Konstanten, 
so findet man durch Einsetzen in (4), daß (60a) verallgemeinert 
werden kann in: 


‘ist eine beliebige Konstante. 


1) Ann. d. Phys. [5] 37. S. 98. 1940; GI. (26). 
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Ein Integral, welches eine willkürliche Funktion enthält, läßt 
sich aus (60b) dadurch gewinnen, daß man r als Funktion von ¢ 
ansetzt und die Einhüllende der Flächenschar 7 = n(z, t, &) aufsucht}, 
Man hat dann neben (60b) die Beziehungen: 


et+e\ _ 1 ate dt 
(61) = T(é); 0= (FS Ge); «const, 


Für jede Wahl der Funktion ze) werden e und z wegen = m- — 


durch z,t ausdrückbar. Wird r(e) willkürlich gelassen, so kommt 
daher durch r eine willkürliche Funktion von 2, t in (60b) hinein, 
Daß (60b) auch unter den Bedingungen (61), also bei nicht- 
konstantem ¢ und r, noch eine Lösung von (4) ist, kann durch 
Einsetzen in (4) bestätigt werden. Wegen der zwei Vorzeichen- 
möglichkeiten von £ kann man für das obere Vorzeichen eine andere 
willkürliche Funktion r(e) wählen als für das untere; man bekommt 
so zwei Lösungen aus (60b), von denen jede mit einer willkürlichen 
Funktion ausgestattet ist. 

Das allgemeine Integral hat man damit nicht; denn die zwei 
eben genannten Lösungen können nicht additiv zu einer zusammen- 
gesetzt werden, die zwei willkürliche Funktionen enthielte; die 
Differentialgleichungen (4) sind ja nicht linear. Es handelt sich 
vielmehr um die allgemeinste Lösung vom Typ &=f(n), also um 
das Integral?) 

t(n + 1) n + 1) 


(62) q F (2 + 2Vn i— nt ) é= + (7 (n- 


e(n+1) 


wie leicht zu beweisen ist und wie schon der Zusammenhang (60b) 
zwischen £ und 7 vermuten läßt. Aus (60b) und (61) folgt 
nämlich, daß: 
(m+) (n + 1) 
+ —— 


——- 
(62a) | 2Yn -1 2Vn-1 2Vn-1 


| 


eine willkürliche Funktion g(e) von & allein ist. Daher ist ¢ eine 
willkiirliche Funktion von ». Nach (62a) ist aber auch: 


(n— 


t(e\(n +1), 
- (7 = gle) 
2yn-ı N gie), 
1) Die Betrachtungen dieses Absatzes entnehme ich einem Brief von 
Herrn Sauter an mich, mit seiner freundlichen Einwilligung. 
2) Es ist in Ann. d. Phys. [5} 37. S$. 98, Gl. (28) angegeben; dort steht 
ı(n +1) 


eine Konstante A an Stelle von i 
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also auch 7 eine willkürliche Funktion von e, mithin auch von a: 
n=Fio»). & bestimmt sich aus (60b): =+ ( _ Damit 


ist der Beweis fertig. 

Daß das Verfahren des vorletzten Absatzes, die Einhüllende einer 
einparametrigen Integralflächenschar aufzusuchen, hier zum Ziel führt 
ist nicht selbstverständlich, ist vielmehr ein seltener Gliicksfall. Denn 
die Bedingungen, [hier (61)], denen der Parameter (hier e) zu genügen 
hat, damit die 7-Fläche Einhüllende der Integralflächenschar (von 
konstantem &) wird, sind im allgemeinen nicht zugleich mit den Be- 
dingungen erfüllbar, die erfüllt sein müssen, damit die durch r=rie) 
entstehende &-Fliche Integralfliche wird. 

In dem oben betrachteten Fall aber sind die Verhältnisse sehr 
einfach: Durch rt = r(e) werden aus den 7- und &£-Flächen einpara- 
metrige Flächenscharen (mit dem Parameter &); die Bedingung 

0 (4 + -0 

Oe\tit 
ist die Bedingung für die Einhüllende der einparametrigen „-Flächen- 
schar, aber zugleich die Bedingung für die Einhüllende der einpara- 
metrigen £-Flächenschar, weil nach (60b) &=+ 4 + const gilt. (Die 
Einhüllenden sind ihrer Natur nach selbst wieder Integraltlächen). 

Bei den Zylinder- und Kugelwellen liegt das Problem anders, 
aus zwei Gründen. Erstens ist die singuläre Lösung (22a) nicht 
verwendbar, weil sie nur eine Konstante enthält (r); die Erweiterung 
r—>r+ const ist hier natürlich nicht möglich; die Erweiterung 
auf die Form £=+ » + const aber auch nicht, wie das Einsetzen 
in (4) zeigt. Der Ansatz &=+n-+ const führt übrigens zwangs- 
läufig auf (22a). Zweitens geben die Lösungen (17) und (52) den 
Zusammenhang: = m7; es ist daher nicht zu erwarten, daß die 
Einhüllenden einparametriger Scharen von 7- und &-Flächen durch 
dieselbe Bedingung geliefert werden. 

Durch geeignete Verallgemeinerung erhält man aber grund- 
sätzlich die Möglichkeit, aus dem vollständigen Integral das allgemeine 
zu finden; leider besteht diese Möglichkeit nur theoretisch, denn sie 
führt auf die Integration eines Systems von zwei simultanen partiellen 
Differentialgleichungen erster Ordnung, die in unserem Fall schwieriger 
sein dürfte als die Lösung von (4). Der Weg zu diesen Gleichungen 
wird im Anhang gezeigt; die Methode ist angegeben im Lehrbuch 
von A.R. Forsyth'. Für k = 1 gibt die Anwendung der Methode 
auf das Integral (60b) die Gl. (61) und damit das Integral (62), wie 
wir am Schluß des Anhangs zeigen werden. 


1) A. R. Forsyth, Theory of Differential Equations, Part IV, Vol. V, 
Cambridge, Univ. Press 1906, S. 419. 
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Anhang 


Über die Ableitung des allgemeinen Integrals 
aus dem vollständigen Integral (52) 


Die Gl. (52) schreiben wir um in’): 


J * de 
0% B, ’ 
913 

die Integrationsvariable ist in dem Integral in (a) x genannt; q enthält 
die Integrationskonstante B von (51) und außerdem den Parameter «, 
wie man an (51) sieht; k und n sind bei gegebenem physikalischem 
Problem feste Zahlen. In der Funktion ¢, die wir in (c) eingeführt 
haben, kommen g, q und « explizit vor, und q enthält selbst wieder g, 
B und «. 

Dies ist für das Folgende wichtig, denn die Größen B, C, r, a, 
die in der bisherigen Behandlung als konstant galten, sollen jetzt 
als Funktionen von r,t — werden. Die Methode besteht 
darin: Damit die Gl. (4) durch (a), (b), (c) erfüllt werden, sollen die 
Zusatzglieder verschwinden, die zu £, E., n,, n, bei der Differentiation 
nacht und r dadurch hinzukommen, daß jetzt die „Konstanten“ B, C,t,« 
von r,t abhängen. Das gibt partielle Differentialgleichungen fiir B,C, r,«. 

Aus (a) folgt durch Differentiation: 


d dC d 
(d) „tier da + — —J,dB — J,de, 


wobei wir die Abkürzungen eingeführt haben: 

oq d 2: 
aus (c) folgt durch en Differentiation: 


dr dt +dı , 
t+r (38 + dp 


Olgs (da j 
| + (24 da) + de; 
(b) gibt: 

(g) 


(f) 


1) Mit der Benutzung von (52) ist die Voraussetzung « + — 1 verbunden. 
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dg kann aus (d) entnommen und in (f), (g) eingesetzt werden. Die 
übrigen Differentiale denken wir durch dr und dt dargestellt und 
erhalten so aus (f) die zwei Gleichungen: 


t+r + (538 04 


N C, 
B, + J,«,| 


Olgt Olgt 

(33 B,+ $4 «,) + da 
ölg: 

+ ( oq 4-9 


(g) liefert ebenso die zwei Gleichungen: 


&= +J, B,+J,a, + $1, : 


Das Gleichungssystem (4) ist jes dann erfüllt, wenn wir 


in (f’)—(g”) die Glieder zum Verschwinden bringen, die von der 
Differentiation der „Konstanten“ herrühren; dabei ist zu beachten, 


daß der Faktor es: + lg: 94 formal genau dasselbe Aussehen 


hat, wie wenn er für das vollständige Integral berechnet würde, also 
mit konstantem B,C,r,«; denn in dem genannten Ausdruck wird 
bei konstantem B, « differentiiert. Wir könnten ihn mit 


bezeichnen; sein Wert kann aus (50d) entnommen werden und ist vor 
(ö3c) angegeben. — Da die so verkürzten GI. (f’)—(g”) die Differential- 
gleichungen (4) erfüllen, wenn die B, C, r, « konstant gelassen werden 
(rgl. § 6), so tun sie es auch bei veranderlichem B,C,r,« noch. 
Aus (f’), (f”) schließen wir so auf die Bedingungen: 
ölg: 


0=— + + aq #) 


| 
| a, Igr - +5, B,, + 


Olgt (dq 07 
(35 2.. + 


| 
| + 
Nr 1 
(f”) | 1 rtd, «,| 
ölg: ($3 
| + + 
| N 
dlg: dig (z) 
di =— ag 
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die Bezeichnüng r, , soll bedeuten, daß die Gleichung sowohl für 
die Differentialquotienten nach t, wie auch für die nach r richtig 
ist; die Gleichungen tragen die Nummern (h,) (für Differentiation 
nach ¢), und (h,) (für Differentiation nach r). (g’), (g”) geben die 
Bedingungen: 


C +d, 


0=a,,lgr— 
Wegen (i, ,) fallen die eckigen Klammern in (h, ,) weg. Die so 
übrigbleibenden 4 Bedingungen können auf verschiedene Weise 
erfüllt sein. 

1. Wir nehmen an, daß « und r voneinander unabhängig seien, 
dann können wir sie als unabhängige Veränderliche an Stelle von r,t 
wählen. Aus den vorigen 4 Bedingungen wird so: 


oq 


dq , Alg: ) 
| 


r 


ö 
k (B, Thr + B, ) 
9! 


C 
O= a, Igr—— 


J, (B, B,«, + J, 
Die triviale Möglichkeit ¢,,= 0, t,,= 0, die wieder auf das voll- 
ständige Integral zurückführt, interessiert uns nicht. In (k, ,) haben 
wir dann zwei lineare Gleichungen für die Faktoren von r, ,, «, ,, 
die in beiden Gleichungen ja dieselben sind. Es gibt die zwei Möglich- 
keiten: Entweder sind diese Faktoren gleich Null; oder es ist die 
Determinante ihrer Koeffizienten gleich Null, dies ist aber die 
ö(t, a) 
d(r, t) 
gegen der Voraussetzung 1. Eine ähnliche Betrachtung läßt sich 
auf (l,,,) anwenden. So kommt im ganzen als erste Möglichkeit: 


dq oB B; (5% B.+ 


Funktionaldeterminante ; ist sie Null, so wird 7 = r(@), ent- 


C, 
|o=Igr—% +J,B,+J.; 0= +J,B.. 


T 


Nach (a) und (c) ist m als Funktion von 1, t,r, B, C, @ grundsätzlich 
also grundsätzlich durch diese Veränderlichen ausdrückbar. Mit 
Hilfe von zweien der Gl.(I) kann man daher r und ¢ einzeln als 
Funktionen von B,C,«,r und von den Differentialquotienten von B,C 
nach @ und r darstellen. Einsetzen dieser Ausdrücke in die übrigen 


berechenbar; die Faktoren 


K. 
Zwe 
glei 
Ihr 
ken 
als. 
in £ 
sich 
gen 
sche 
auc 
gefi 
Var 
Vor 
abhi 
Doc] 
und 
(m,, 
{Dh 
Nac] 

(Il) 
Wen 
awe wihl 
gi. 
| 
| 
Bezic 
durel 
eine 
die ] 
B, B 


K. Bechert. Ausbreitung von Zylinder- und Kugelwellen usw. 199 


wei Gleichungen gibt dann zwei simultane partielle Differential- 
gleichungen erster Ordnung fiir B, C als Funktion von @ und r. 
Ihre Lösung enthält zwei willkürliche Funktionen von «, r; man 
kennt damit B, C in Abhängigkeit von a, r und kann r und t einzeln 
als Funktionen von «, t darstellen. So sind auch «, r und damit B, C 
in Abhängigkeit von r,t bekannt und die Aufgabe ist gelöst. Wenn 
sich das Programm praktisch durchführen ließe, wäre damit das all- 
gemeine Integral gefunden. Daß es nicht durchführbar ist, sieht man 
schon an den Gl. (a)—(c). — Bemerkenswert scheint mir, daß man 
auch hier wieder auf die „Umkehrung“ (wie beim ebenen Problem) 
geführt wird: man würde zunächst r, t in Abhängigkeit von anderen 
Variablen bekommen. 

2. Die andere Möglichkeit: + = r(«) fällt nicht mehr unter die 
Voraussetzung 1. Wenn r=r(«), kann man «& und r nicht als un- 
abhängige Veränderliche benützen; die G1. (h), (i) gelten aber natürlich 
noch, Wir wählen jetzt @ und C als unabhängige Veränderliche 
und erhalten aus (h), (i): 

Oleg: 
+73 


ög ölg: 


0=— 


(B, + Be C,,r) 


+ I, (By tyr + + Jy 


O=e,,lgr— 


Nach demselben Verfahren wie in 1. schlieBen wir auf: 


ölg: Olg:. 


+d, Be. 


Wenn 3, = U ist, so gilt entweder: B ist konstant und fest ge- 


wählt; dann ist aber auch B. = 0, was einen Widerspruch zur letzten 
Gl. (II) in sich schließt. Oder: es Egg eine Gleichung zwischen B, 


4, @ allein. Wegen (a), (c) und = 0 kann g aber ausgedrückt 


werden durch r, t, C, «, r. Daher haben wir eine Beziehung 
zwischen B, r, t, C, «, r; daraus kann r als Funktion von t, C,«@,r, B 
berechnet werden. Einsetzen in die dritte Gleichung etwa gibt eine 
Beziehung zwischen t, C, a, t, B, B,. Man kann demnach t ausdrücken 
durch ©, «@,r, B, B,. Einsetzen von r,t in die erste Gleichung gibt 
eine Beziehung (A) “zwischen C, a, t, B, B,, t,; Einsetzen von r,t in 
die letzte Gleichung liefert einen Zusammenhang zwischen C, «, r, 
B, B,, Be. Integration dieser Gleichung gibt B als Funktion von « 
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und ©, wobei noch eine Funktion von «, außerdem C und ze) will. 
kürlich bleiben. Das so gefundene B wird aber im allgemeinen die 
vorhergehende Gl. (A) nicht erfüllen; aus dieser müßte ja r be. 
rechenbar sein und dürfte nur von & abhängen, was eine Ein. 
schränkung für B bedeutet. Sollte diese einschränkende Bedingung 
identisch erfüllt sein, so würde sich ein Integral von (4) ergeben, das 
eine willkürliche Funktion von r, t enthielte. — Wenn oe 
erhält man einen Zusammenhang zwischen g und a, also zwischen r,t 
B, C, a, t. Die eben gemachten Schlüsse gelten ungeändert auch fiir 
diesen Fall. Da die wirkliche Durchführung der Rechnungen wohl 
unmöglich ist, wird man so nicht entscheiden können, ob (II) mögliche 
Lösungen gibt. — Die Gl. (m), (n) lassen — ähnlich wie (k), (I) in 1.— 
die Möglichkeit offen, daß neben rt = r(«) auch noch 
3.C=(C(e) gilt. Dann folgt aus (h), (i): 

öq 


= 0 ist, 


(04,2) 


Das gibt: 


ögq Oa 
+J,; O=J,. 


| 0= lgr _ 
Durch B = const wird die zweite und die vierte Gleichung identisch 
erfüllt, vgl. (e). Aus der ersten Gleichung ist dann r als Funktion 
von t, C, a, t, t, berechenbar; Einsetzen in die dritte Gleichung gibt 
eine Beziehung zwischen t,C, «,r,r,,C,. Hier muß entweder ¢ ganz 
herausfallen, was eine Bedingung ist, die im allgemeinen nicht erfüllt 
sein wird, oder es wird « Funktion von t allein; dann muß r au 
der ersten Gleichung herausfallen, was wiederum im allgemeinen 
nicht als zutreffend erwartet werden kann’), Sollte eine der genannten 
Bedingungen erfüllt sein, so würde man im ersten Fall etwa C(@ 
willkürlich wählen können, daraus z(«) finden und schließlich « in 
Abhängigkeit von r,t bekommen. Im zweiten Fall wäre t =t(@ 
willkürlich wählbar und daraus C(«), r(«) zu berechnen. In beiden 
Fällen würde sich ein Integral von (4) mit einer willkürlichen Funktion 
ergeben. — Zur Erfüllung von (o), (p) hat man aber noch die Mög- 
lichkeit: 


1) Natürlich ist auch der Fall denkbar, wo r und ¢ ihre Rolle vertauscht 
haben. 
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4. B= B(e); zugleich gilt gemäß den gemachten Voraus- 
setzungen: r=r(e)}; C= (Cie). Die Gl. (hi, (i) geben jetzt: 


| O=lgr— 45, B, +4,; 


(IV) 


hier kann derselbe Weg gegangen werden wie in 3. Berechnen wir 
nämlich r aus der ersten Gleichung als Funktion von t,C,a@,t, B, B.,t,, 
und setzen in die zweite Gleichung ein, so muß wieder gelten: ent- 
weder fällt ¢ heraus oder es wird @ Funktion von t allein, dann 
muß dafür r aus der ersten Gleichung herausfallen. Im allgemeinen 
wird keine dieser Bedingungen erfüllt sein. 

Bis jetzt haben wir a +— 1 vorausgesetzt. « =— 1 verlangt 
die Verwendung von (55), (56); wir haben nur 3 Konstanten zur Ver- 
fügung; also kein vollständiges Integral, das hindert die Anwendung 
der Methode aber nicht. Man kann leicht die Beziehungen ableiten, 
die den Gl. (h), (i) entsprechen. Es gibt jetzt nur 3 Fälle der Erfüll- 
barkeit dieser Bedingungen, und da wir nur 3 Konstanten haben, so 
sind diese 3 Fälle im allgemeinen nicht als realisiert zu erwarten. 
Das kann man sich genau so überlegen wie in den vorhergehenden 
Fällen 2.—4. bei e+-— 1. 

Das Gesamtergebnis ist also: Der Weg 1. (« + — 1) ist grund- 
sätzlich gangbar, und liefert das allgemeine Integral; in unserem Fall 
verlangt das Programm die Integration eines Differentialgleichungs- 
systems, das vermutlich komplizierter ist als das ursprünglich ge- 
gebene System (4). Die übrigen Möglichkeiten 2.—4. geben im all- 
gemeinen keine Lösungen, die über das vollständige Integral hinaus- 
gingen. 

Einen solchen Ausnahmefall, wo eine der zuletzt genannten 
Möglichkeiten ein brauchbares Ergebnis liefert, wollen wir zum 
Schluß behandeln. Wir werden zeigen, daß bei den ebenen Wellen 
die Anwendung der Methode auf das Integral (60b) notwendig auf (61) 
führt. (60b) ist kein vollständiges Integral, es kommen nur drei 
willkürliche Konstanten darin vor; wir bekommen auch nicht das 
allgemeine Integral, aber eines mit einer willkürlichen Funktion. 

Man hat nach (60b): 


2Vr-1 & 


2Vn—ljlt+e, (x+8)1, 
{ - 14 2 


§=+ (n ; 


n-1 
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Wir verlangen daher, daß: 


2Vn-1 &, (vw + 2 
(41,2) n+1 (¢ + 1)? 0: 


2) 


von (4,,,) bleibt so nur: 
(e+E&T,x Pr 
(8), 2) 


Da die Koeffizienten von &,, und z, , in (8, ,) nicht alle Null sein 


können — Fall 1. scheidet also aus —, so bleibt nur r = r(9, 
was 2. entspricht. Aus den Gl. (s) wird so: 


=(; 
€ 


dazu kommt noch r = r(e) und (r, ,), also « = const; das ist genau 
der Inhalt der Gl. (61). 


GieBen, Institut fiir theoretische Physik, 23. Januar 1941. 


(Eingegangen 26. Januar 1941) 
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Dielektrische Nachwirkungserscheinungen 
an Phosphoren') 


Von Hans Ruffler 
(Mit 4 Abbildungen) 


Die Untersuchung der elektrischen Eigenschaften phosphoreszenz- 
fähiger Stoffe ist ein frühzeitig insbesondere von Lenard erkannter 
Weg, zur Kenntnis des Mechanismus des Phosphoreszenzvorganges 
und der Molekülkomplexe vorzuschreiten, die als Leuchtzentren für 
das Erscheinungsgebiet maßgebend sind. Faßt man diese Zentren 
nach Lenard?) als elektrische Oszillatoren und damit als elektrische 
Dipole auf, so muß deren Anzahl, elektrisches Moment und Beweg- 
lichkeit das dielektrische Verhalten dieser Stoffe bestimmen und 
andererseits aus letzterem erkennbar werden. Die in dieser Richtung 
geführten Untersuchungen, wie sie auf Dielektrika allgemein an- 
wendbar sind’), erscheinen bei den Phosphoren deshalb besonders 
aussichtsreich, weil man es hier in der Hand hat, die verantwortlichen 
Molekülkomplexe in einfacher und eindeutiger Weise nach Zahl und 
Beschaffenheit durch ein erprobtes Präparationsverfahren herzustellen. 

Die vorliegende Arbeit befaßt sich mit den dielektrischen 
Nachwirkungserscheinungen, wie sie bei statischen Messungen im 
zeitlichen Gang der Dielektrizitätskonstanten (DK.) bei CdS—Cu-Phos- 
phoren verschiedenen Kupfergehalts kenntlich werden. Ihre Ergebnisse 
entsprechen der Vorstellung vom Vorhandensein von Dipolen, deren 
Beweglichkeit und elektrisches Moment in definierter Weise von 
der Struktur und dem Erregungszustand des Phosphors abhängt. 


Meßmethode 


Die Beobachtung beruht auf der kurzzeitigen Verfolgung des 
Nachladestroms in einem mit dem Phosphor beschickten Konden- 
sator nach Anlegen einer Spannung. Das Verfahren ist im wesent- 
lichen aus älteren Untersuchungen des Verhaltens unvollkommener 


1) D16. 

2) P. Lenard, Verh. des Naturhist. Med. Vereins Heidelberg 10. S. 7 
1909; Handb. der Exp.-Physik S. 98. 1928. 

3) Vgl.E.v.Schweidler, Ann. d. Phys. 24. 8. 711. 1907; K. W. Wagner, 
Ann. d. Phys. 40. S. 817. 1913. 
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Dielektrika bekannt'). Deshalb genügen hier einige kurze, den vor- 
liegenden Besonderheiten gewidmete Hinweise. Die MeBanordnung 
zeigt schematisch die Abb. 1. Der Kondensator C enthält den zu 
untersuchenden Phosphor; er ist so gebaut, daß er ohne Füllung keine 
dielektrische Nachwirkung zeigt, 
Seine Leerkapazität beträgt etwa 
20 pF. Da der Ausschluß von 
Feuchtigkeit unerläßlich ist, ist er 
in ein mit Phosphorpentoxyd be- 
schicktes evakuierbares Gefäß ein- 
gebaut, durch dessen Glaswand hin- 
Abb. 1. Schaltbild durch er belichtet werden kann. Als 
der MeBanordnung Meßinstrument wird ein Quadrant- 
elektrometer benutzt, dem eine 
Ballastkapazität K von etwa 0,1 „F parallel geschaltet ist. Der 
Schalter S,’, der die Kapazität A kurzzuschalten gestattet, ist gleich- 
zeitig der Erdschalter des Elektrometers, der vor jeder Messung zu 
öffnen ist. 

Der zeitliche Gang der Aufladung von C wird nun in folgender 
Weise gemessen: S,’ ist dauernd geöffnet und außer Benutzung. 
Schalter U liegt am Widerstand R, S, und S, sind geschlossen; 
eine Aufladung von C erfolgt noch nicht. Wird S, geöffnet, so 
fließt der jetzt einsetzende Ladestrom von C über S, nach K so lange, 
bis S, geöffnet wird. Aus der sich einstellenden Spannung am 
Elektrometer findet man mittels der Kapazität von K die auf- 
genommene Ladung, der eine durchschnittliche Dielektrizitätskon- 
stante des Phosphors zugeordnet werden kann. Die Variation der 
Aufladezeit erfolgt in kleinen Intervallen von Zehntel-Sekunden mit 
Hilfe des ballistischen Pendels, das nacheinander die Öffnung vonS, 
und S, betätigt. 

Da die benutzten Phosphore neben der dielektrischen Polarisation 
eine merkliche Leitfähigkeit besitzen, so muß diese getrennt gemessen 
und vom Gesamteffekt in Abzug gebracht werden. Für diese Messung 
wird der Schalter S,’ benutzt, während S, dauernd geöffnet außer 
Funktion bleibt. Man legt U an R und schließt S,’ und S;, so daß der 
Kondensator C gegen Erde sich auflid. Nach gemessener Zeit 
wird wieder mittels des Pendels S,’ und kurz darnach S, geöffnet, 
so daß im Elektrometer diejenige Ladung angezeigt wird, die jeweils 
in der Zeitspanne zwischen dem Öffnen von S, und S, noch über C 


1) Vgl. F. Tank, Ann. d. Phys. 48. S. 307. 1915; O. Gayling von Alt- 
heim, Ann. d. Phys. [5] 35. $. 417. 1939. 


get 
hal 
der 
die 
die 
wal 
nat 
30( 
der 
sin 
Pri 
Cu- 
Phe 
din; 
der 
Err 
diel 
ans 
szis 
jew: 
Phe 
als 
der 
ten 
che 
gem 
fach 
Bez 
der) 
mög 
0,01 
geri 
ansı 
ursa 
bem 
aus; 


H. Ruffler. Dielektrische Nachwirkungserscheinungen usw. 205 


geflossen ist. Man beobachtet, daß diese Ladungsmenge bei Konstant- 
haltung der genannten Zeitspanne etwa 10 Minuten nach Anlegen 
der Spannung E an C konstant wird und daher dann als Maß für 
die Ohmsche Leitung im Phosphor angesehen werden kann. Falls 
dieser Strom nach Art eines Sekundärstroms zeitlich inkonstant 
wäre, was hier nicht untersucht werden konnte, wäre das Elimi- 
nationsverfahren nur in Näherung zutreffend. 

Die zur Messung benutzten Feldstärken werden von 20 bis 
3000 Volt/em im Phosphorkondensator variiert. Es zeigt sich, daß 
deren Höhe den Gang der e-Werte nicht beeinflußt. 


Meßergebnisse 


Alle Ergebnisse beziehen sich auf CdS—Cu-Phosphore, die 
sämtlich aus dem gleichen Grundmaterial unter genau denselben 
Präparationsbedingungen gewonnen sind. Variiert ist lediglich der 
Uu-Gehalt in den weiten Grenzen von 0,01—0,00001 g Cu im Gramm 
Phosphor. Zum Vergleich kommt noch das unter den gleichen Be- 
dingungen geglühte reine Grundmaterial CdS. Untersucht wird sowohl 
der Dunkeleftekt am unerregten Phosphor als auch der Einfluß der 
Erregung. 

1. Einfluß der Cu-Konzentration. — Der zeitliche Gang der 
dielektrischen Nachwirkung wird durch Abb. 2 in der Weise ver- 
anschaulicht, daß als Ab- 
szisse die Aufladedauer des 9 Up 
jeweils zuvor entladenen 
Phosphorkondensators und 
als Ordinate ein Relativwert 
der Dielektrizitätskonstan- 
ten & verzeichnet wird, wel- 
cher der am Elektrometer 


gemessenen Aufladung ein- B BH Z 


fach ti ist. 
proportional ist. Als Abb. 2. Einfluß der Cu-Konzentration 


Bezugswert der DK. wird auf die dielektrische Nachwirkung 
derjenige &, gewählt, der 


möglichst unmittelbar nach der Spannungsanlegung, praktisch nach 
0,01 Sek. beobachtet werden kann. 

Wie man erkennt, ist die Nachwirkung beim Grundmaterial sehr 
gering. Ihre erhebliche Steigerung durch Kupferzusatz deutet sehr 
anschaulich auf das durch die Ausbildung der Leuchtzentren ver- 
ursachte erhebliche Anwachsen der Zahl trägerer Dipole hin. Sehr 
bemerkenswert ist die Feststellung, daß die Nachwirkung um so 
ausgeprägter ist, je geringer die Metallkonzentration, und daß die, 
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größte benutzte Cu-Menge die geringsten Veränderungen gegenüber 
dem reinen Grundmaterial hervorruft. Dieses Verhalten ist verständ- 
lich und ein Hinweis auf die schon von Lenard festgestellte Tat- 
sache, daß mit steigendem Metallgehalt die Zentrengröße eines Phos- 
phors abnimmt, so daß offenbar ihre Orientierungsmöglichkeit im 
elektrischen Feld erleichtert wird. Da alle mitgeteilten Kurven auf 
gleichen Ausgangswert &, bezogen sind, machen sie naturgemäß nur 
eine Aussage über das Zurücktreten von Dipolen großer Zeitkonstanten 
(im Sinne K. W. Wagners) mit wachsendem Metallgehalt. Dagegen 
bleiben die Gesamtzahl der Dipole und ihre elektrischen Momente 
unbekannt, solange die Kenntnis der in Betracht kommenden Absolut- 
werte von s, deren einwandfreie Ermittlung nicht ohne Schwierig- 
keiten ist, aussteht. Auch wäre es wohl erforderlich, die Nachwirkung 
zeitlich noch erheblich weiter zu verfolgen. 
2. Einfluß der Temperatur. — Wie bei den Nachwirkungsvor- 
gängen anderer Dielektrika findet man auch bei den Phosphoren 
eine Temperaturabhängig- 
keit. In Abb. 3 findet man 
u u al die Meßergebnisse für die 
y 4 beiden Phosphore mit den 
Va extremen Cu-Gehalten von 
0,01 und 0,00001 g im 
ww. Gramm Phosphor und für 
ai 3 Temperaturen zusam- 
a mengefaBt. Die Nach- 
ie wirkung ist bei allen 
Temperaturen um so grö- 
Abb. 3. Einfluß der Temperatur Ber, je geringer die Metall- 


auf die dielektrische Nachwirkung konzentration. Sie nimmt 
außerdem mit wachsender 


Temperatur stark zu und deutet insbesondere eine vermehrte Mit- 
wirkung großer Dipole an — Abnahme ihrer Zeitkonstanten — 
die bei gewöhnlicher Temperatur ihrer geringen Beweglichkeit wegen 
innerhalb der verfolgten Nachwirkdauern nur einen sehr geringen 
Beitrag liefern. Der weitere Aufstieg der Kurven konnte auch noch 
bei einem versuchsweisen Erhitzen einiger Phosphore auf etwa 200°C 
beobachtet werden; ein Sättigungswert von & war also auch hier 
noch nicht angezeigt. 

3. Einfluß der Belichtung. — Die bei ZnS-Phosphoren lange 
bekannte Steigerung der DK. durch Belichtung tritt auch bei den 
CdS-Phosphoren in ausgeprägtem Maße auf. Die vorliegenden Be- 
obachtungen zeigen, daß auch die Nachwirkungserscheinung durch 
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Belichtung eine Steigerung erfährt. Die Erregung, über deren Er- 
gebnis hier berichtet werden soll, erfolgt mit weißem Licht, dessen 
roter und ultraroter Anteil durch eine zwischengeschaltete Kupfer- 
sulfatlösung absorbiert wird, um die langwellige Tilgung auszuschalten. 
Gemessen wird die in bestimmter variierter Zeit über den Phosphor- 
kondensator geflossene Elektrizitätsmenge im unbelichteten und im 
belichteten Zustand. Um den sich ergebenden Unterschied klar 
hervortreten zu lassen, sollen in Abb. 4 nicht die den Verlauf dar- 
stellenden Einzelkurven, sondern deren Differenzen, die also den 
reinen Belichtungseinfluß ver- 
anschaulichen, verzeichnet | 
werden. Die Ordinatenwerte | | 
der Kurven geben, im loga- 
rithmischen MaBstab aufge- 
tragen, denjenigen zusätz- 
lichen Teil der Kondensator- 
aufladung in Coulomb, wel- 
cher ausschließlich als Folge 
der Belichtung in den auf 
der Abszisse genannten Zeiten 
auftritt. Wir beschränken 
uns auch hier auf die beiden oe | 


extremen Phosphorproben: a 


Abb. 4. Einfluß der Belichtung 
— auf die dielektrische Nachwirkung 


liegenden Metallgehalten zu- bei verschiedenen Temperaturen 
gehörigen Kurven würden 
im Zwischengebiet zwischen den mitgeteilten liegen. 

Man erkennt, daß der Einfluß der Belichtung auf die dielektrische 
Nachwirkung bei hohem Metallgehalt des Phosphors erheblich geringer 
ist als bei geringem Gehalt. Diese Abhängigkeit ist in Überein- 
stimmung mit derjenigen der optischen Erregbarkeit der Phosphore; 
sie weist also ebenfalls auf die Beteiligung der Leuchtzentren am 
Effekt hin. Beachtenswert ist die Beobachtung, daß der Temperatur- 
einfluß beim belichteten Phosphor!) erheblich größer ist als beim 
unbelichteten (Abb. 3), und daß er bei hoher Metallkonzentration 
stärker hervortritt, während dies beim Dunkeleffekt nicht der Fall 
ist, Es werden damit Vorgänge angedeutet, deren eingehendere 
Verfolgung unter Variation der Erregungsweise im Zusammenhang 
mit der Verfolgung der jeweiligen optischen Eigenschaften geeignet 


hung 
S_ 


& 


1) Uber den Temperaturgang der im Hochfrequenzfeld gemessenen DK. 
eines ZnS-Cu-Phosphors vgl. R. Wilde, Ztschr. f. Phys. 15. 8. 350. 1923. 
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erscheint, einer vertieften Kenntnis des Zentrenaufbaues der Phos- 
phore näher zu kommen. Die weiterhin zu erstrebende absolute 
quantitative Verwertung der Meßergebnisse macht schließlich die 
Bestimmung der wirksamen Stoffmenge und optischen Energien er- 
forderlich. 

Zusammenfassung 


Es wird die dielektrische Nachwirkung an CdS—Cu-Phosphoren 
variierten Cu-Gehalts im unerregten und erregten Zustand bei ver- 
schiedenen Temperaturen untersucht. Dabei findet sich: 

1. Bei konstanter Temperatur nimmt die Größe der Nachwirkung 
mit abnehmender Metallkonzentration zu. 

2. Die Belichtung führt zu einer erheblichen Steigerung des 
Nachladestroms. Diese Steigerung: ist bei geringer Metallkonzen- 
tration größer als bei hoher Konzentration. 

3. Temperatursteigerung vergrößert die Nachwirkung. Diese 
Vergrößerung ist beim unerregten Phosphor nicht wesentlich vom 
Metallgehalt abhängig. Durch die Belichtung wird der Temperatur- 
einfluß auf den Nachladestrom erheblich gesteigert, und zwar um 
so mehr, je höher die Metallkonzentration ‘ist. 

Führt man die typischen dielektrischen Eigenschaften der Phos- 
phore auf das Vorhandensein elektrischer Dipole zurück, die als 
Leuchtzentren gleichzeitig das optische Verhalten derselben bestimmen, 
so erscheinen Beobachtungen der vorstehenden Art geeignet, der 
Kenntnis der Zentrenstruktur der Phosphore näher zu kommen. 
Die dabei zu erstrebenden mehr quantitativen Festlegungen werden 
allerdings vor allem zu entscheiden haben, wie weit die Dipolvor- 
stellung allein genügen kann, um dem Gesamt-Verhalten dieser in- 
homogenen Dielektrika ausreichend gerecht zu werden. 


Herrn Prof. Dr. A. Becker, dem Direktor des Institutes, sowie 
Herrn Prof. Dr. L. Wesch danke ich herzlichst für die stete Hilfe 
und Förderung, die sie dieser Arbeit angedeihen ließen. Ebenso 
herzlich danke ich meinem Freunde Dr.K. Kamm, der mir in freund- 
schaftlicher Weise die hier untersuchten Phosphore präparierte. 


Heidelberg, Philipp Lenard-Institut der Universität. 


(Eingegangen 21. Januar 1941) 
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Schallabsorption in Mischungen einatomiger Gase 
Von Max Kohler 


Die Schallabsorption ist in reinen einatomigen Gasen nur durch 
die innere Reibung und die Wärmeleitung bedingt. Die grundlegende 
theoretische Untersuchung über diesen Gegenstand stammt von 
Kirchhoff!) aus dem Jahre 1868. Für reine einatomige Gase stellen 
die Kirchhoffschen Formeln den theoretischen Sachverhalt er- 
schöpfend dar. Die experimentellen Bestimmungen der Schall- 
absorption an Kdelgasen bestätigen zwar die aus Kirchhoffs 
Formeln folgende Druck- und Frequenzabhängigkeit, ergaben aber 
etwas zu hohe Zahlwerte. Ein Versuch zur Erklärung dieser Dis- 
krepanz wurde von K. F. Herzfeld?) unternommen, indem er bei 
der Ableitung der hydrodynamischen Grundgleichungen im Rahmen 
der kinetischen Gastheorie die Annäherung einen Schritt weiter treibt. 
Die Frage nach der Schallabsorption in Gasgemischen wird deshalb 
nicht durch Kirchhoffs berühmte Abhandlung beantwortet, weil in 
diesem Fall noch Diffusionserscheinungen hinzukommen, die die 
theoretische Behandlung des Problems sehr erschweren. Während 
die durch Reibung und Wärmeleitung auftretenden Schallverluste 
noch den phänomenologischen hydrodynamischen Grundgleichungen 
entnommen werden können, ist dies bei Gasgemischen nicht mehr 
der Fall. Die Bewegungsgleichungen fiireein Gasgemisch können 
nur aus der kinetischen Gastheorie gewonnen werden. Da bei Schall- 
wellen in Gasen Temperaturdifferenzen auftreten, kommt zu der ge- 
wöhnlichen Diffusion noch die Thermodiffusion hinzu. Es ist daher 
zu erwarten, daß der Schallabsorptionskoeffizient in Gasmischungen 
neben dem Reibungskoeffizienten und dem Wärmeleitvermögen auch 
den Diffusions- und Thermodiffusionskoeffizienten enthält. Prinzipiell 
hat man dann die Möglichkeit, durch Schallabsorptionsmessungen an 
Gasmischungen den sonst der Messung nur schwer zugänglichen 
Koeffizienten der Thermodiffusion zu finden. Außerdem erhält man 
Aufschluß über die Schallabsorption in Isotopengemischen. Nun 
liegen theoretische Untersuchungen über Schallabsorption in Gas- 


1) G. Kirchhoff, Pogg. Ann. 134. S. 177. 1868. 
2) K. F. Herzfeld, Ann. d. Phys. [5] 28. S. 465. 1935. 
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gemischen schon seit längerer Zeit von Y. Rocard!) vor. Rocard 
verwendet bei seinen Untersuchungen Bewegungsgleichungen, wie sie 
aus Chapmans?) Arbeiten über die kinetische Theorie der Gas- 
mischungen folgen. Im Falle von Luft findet Rocard, daß die 
Schallabsorption durch Diffusionsvorgänge etwa von derselben Größe 
ist wie die durch Wärmeleitung bedingten Verluste. . Dieses Ergebnis 
beruht zum Teil auf unzulässigen Anwendungen von Chapmans 
Formeln. Wie von Enskog?) gezeigt wurde, sind überdies ver- 
schiedene Ergebnisse Chapmans über Diffusion und Thermodiffusion 
fehlerhaft. Es erschien daher wünschenswert die Fragestellung noch 
einmal aufzugreifen. Zweckmäßig ist es hierbei die Enskogsche‘) 
Form der Gleichungen zu verwenden. Es gelingt dadurch Formeln 
für die Schallabsorption in Gasgemischen zu erhalten, die ohne 
spezialisierende Annahmen über das Wechselwirkungsgesetz zwischen 
den Atomen gewonnen werden. Die neuen Ergebnisse weichen so- 
wohl qualitativ als auch quantitativ in vieler Hinsicht von denen 
Rocards ab. 


$1. Ableitung der Bewegungsgleichungen 


Es bedeute ¢ eine differenzierbare Funktion von z, y, 2, t, v,, vv, 
die als der zahlenmäßige Ausdruck irgendeiner additiven Eigenschaft 
des einzelnen Moleküls gelten kann. Bei der Ableitung der Be- 


wegungsgleichungen gehen wir dann von folgenden Gleichungen aus’): 


dw 
+ ox dy (n, 


gu) 


“te 


Ow 


1) Y. Rocard, Journal de Physique et le Radium 7. S. 426. 1930, ,,Propa- 
gation et absorption du son“ Paris 1935. Hermann et Cie., Act. Scientif. et Industr. 

2) 8S. Chapman, Phil. Trans. Roy. Soc. 217. S. 215. 1917. 

3) D. Enskog, Arkiv för Mat. Astr. och Fysik 16. 8. 57. 1922. 

4) D. Enskog, Diss. Uppsala 1917. 

5) D. Enskog, Diss. S. 13. 
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Hierin ist n, die Zahl von Molekülen der 1. Gasart pro Kubikzenti- 
meter, ”, bedeute die entsprechende Anzahl der 2. Gasart. Die 
Größen v,, v,, ©, sind die Geschwindigkeitskomponenten der Atome. 
Durch Überstreichen ist eine Mittelwertsbildung über die Atome der 
1. Art angedeutet. Die Größen X,, Y,, Z, sind die Komponenten 
der äußeren Kräfte auf die Atome der 1. Art pro Masseneinheit, die 
Größen x, ), 3 die Geschwindigkeitskomponenten v,— u, v, — v, u, 
wobei u, v, w die Komponenten einer zunächst nicht näher definierten 
Geschwindigkeit seien. Außerdem ist: 


co 2x 


ff dbde, 
0 0 


co 2a 


fff do, db de. 
0 0 


Hierin sind f und F die Verteilungsfunktionen fiir das Gas 1 und 
das Gas 2, im übrigen sind die Bezeichnungen so wie bei Enskog 
und Boltzmann?) nach ad Zeit in (la) ist 


eine Abkiirzung fiir t¥ - Eine zu (la)ent- 


sprechende zweite Gleichung läßt sich ? ar die 2. Gasart hinschreiben. 
Sie werde mit (1b) bezeichnet. 

Es sei nun in (la) zunächst: u=u,,v=0,,w = W,, WO U,,%,, %, 
und u,,0,,w, die mittleren Geschwindigkeitskomponenten der Moleküle 
des 1. bzw. 9. Gases sind. Sind m,, m, die Massen der beiden Atom- 
sorten, o, und o, die Partialdichten der beiden Gase, so ist 9, = n, m, 
und 0,=n,m,. In (la) sei außerdem gy = m,. Wir erhalten dann, 
und 4, (m) = 4, m) 


die Kontinuitätsgleichung des 1. Gases. BE ergibt sich die 
Kontinuitätsgleichung des 2: Gases. Sie lautet: 


ö 
(2b) +0 Us + v + St) =0. 


In (2a) und le sind die ses d,/dt und d,/dt Ab- 
kürzungen für: + +u, +, Se +w, und den entsprechenden 


d ou 
9, 1 
(2a) “a +9 (3 de Oy 


Operator, der durch Vertauschung von W, mit W,, 0,, w, 
hervorgeht. 


1) L. a: „Vorlesungen über kinetische Gastheorie“. 
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Nun mögen u,v, w in (la) und (1b) die Komponenten der 
mittleren Geschwindigkeit aller im Volumelement vorhandenen Atome 
sein. Bezeichnen wir diese mit %,, ?,, %,, so ist: Nu, =n, u, +n, u, 
won=n,+Nn, die Gesamtzahl von Atomen pro 
deutet. Entsprechende Gleichungen gelten für v, und w,. Außerdem 
setzen wir in (la) und (1b) g = 1. Es ergeben sich 2 Gleichungen, 
von denen die erste lautet: 


du, On, (u, — Up) 
+n, (2% + ) + 


On, — Wo) 


=e. 


n, (Uy — Vy) 
öy 
Durch Addition der beiden Gleichungen folgt: 


Uo Ow, \ _ 


Ferner nehmen wir nun an, daß u,v, w in (1) die Komponenten der 
Schwerpunktsgeschwindigkeit der in einem Volumelement befindlichen 
Atome sind, die in folgender Weise definiert sind: 9u =o, u, + 0,4, 
mit zwei weiteren analogen Gleichungen für die y- und z-Komponente, 
Außerdem setzen wir in (la)@=m,. Dann folgt: 


09, (u, — 


do, Ou dw u) — v) 
dt +0, (5" + 


öy 
(w, —w) 

Oz 
Für das 2. Gas gilt eine analoge Gleichung. Durch Addition der 
beiden Gleichungen folgt: 


=0. 


do ou Ov , öw\ _ 

In Zukunft bedeuten u,v, w immer die Komponenten der Schwer- 
punktsgeschwindigkeit. Wir setzen: 


(4) 


(1) a2 (1) (1) 


2. 
2 2 


(1) 

Pry = MM, t,9, usw., 
(2) 


(2) 
Pay = N, Mm, 
(1) 


Schließlich sei: p= + . In (la) setzen wir 


m, (v, — u), 


M 
| 
Kir 
sic] 
(da 
Auf 
Sch 
der 
Die 
pun! 
und: 
n,m, po = n,m, 
usw. 
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dann ergibt sich, da 4,, (m, xt) = 0: 


(u—u+ 0; (u, — u) dr by 


zz =. 
ou ou 
= 0, (x, = — 9, (u, — u) 
6 
— 9, — + 4), (m, 


Eine analoge Gleichung gilt fiir das 2. Gas. Durch Addition ergibt 
sich unter Berücksichtigung von 4,,(m,%) + 4, m, = 0. 


d Op: a: 
ft + = 0, X, + 9, X,. 


Dazu kommen noch 2 weitere Gleichungen, die man vollkommen 
analog erhält, und die die Nummern (5b) und (5c) tragen mögen. 
Wir definieren weiter: 


Q, = mm 39) 
Yo” + 3,2) 
n,m, 


Außerdem sei 9 = Q, +Q, und q,=q,” +4q,® usw. Relativ zum 
Schwerpunkt kann bekanntlich die Bewegungsenergie in die Energie 
der geordneten und die der ungeordneten Bewegung zerlegt werden. 
Die Größe Q ist die Letztere pro Volumeinheit. Die Größe q, ist die 
z-Komponente des Wärmeflusses der Mischung relativ zum Schwer- 


punkt. Wir setzen nun in (la) g = 3 m, t® und erhalten: 


| 
dt Ox oy + Oz 


a) du a) dr a Ow (1) ou 


du dw agi! 


+P 
= 0, (u,—%) (x,- 7) +4, 4-9) av) 


+0, (w, — w) (z, - 4, r’). 


e 
(1) 
x: 
z 
n 
ly 
h 
2 
2 0, ’ 
1 
2 
3 9% t,”, 
q 
r 
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Die Größe 4,, (3 r*) verschwindet. In analoger Weise ergibt sich 


fiir das 2. Gas eine entsprechende Gleichung. Da 
) + 4, (4 ) = 0 
ist, ergibt die Addition der beiden Gleichungen: 
dq ou ou Og: dy 
= 0, (4, — WX, + 0, (Uy — WX, + (0%, — Y, + @, (v, — oY, 


+ 9, (w, — w)Z, + 0,(w, — WZ, — + 


Subtrahiert man die mit Q/o multiplizierte Kontinuitätsgleichung (4), 


so 


a( 
+ 


dv Ow 
(u, u) X, + 0, (U, = u)X, + 9, (v, Y, 
+ 04(%, — )Y, — w)Z, + 9, (w, — w)Z,. 
Zu den bisher gewonnenen Gleichungen treten nun noch Beziehungen 
hinzu, die sich nicht aus (la,b) ableiten lassen, sondern die die 
explizite Bestimmung der Verteilungsfunktionen f und F aus da 
Boltzmannschen Gleichungen voraussetzen. Während die Gl. (2) 
bis (6) vollkommen exakt gültig sind, erfordert die Lösung der 
Boltzmannschen Integralgleichungen wesentliche einschränkend 
Voraussetzungen, die den Gültigkeitsbereich der nach den üblichen 
Verfahren gewonnenen Formeln auf das Gebiet der mäßig ver 
dünnten Gase beschränken. Enskog)!) findet in diesem Fall folgende 
Diffusionsgleichungen: 


Op, oT 


mit 2 weiteren Gleichungen für die y- und z-Komponenten, die 
die Nummern (7b) und (7c) tragen sollen. Die Größe x,, hängt 
direkt mit dem gewöhnlichen Diffusionskoeffizienten zusammen: 


(8) D,, =%,06%P- 


1) D. Enskog, Diss. S. 70. 


und s 


Verh: 
atomi 
sonst 
Molel 
Größe 
vorgä 
M= 

gegen 
\ 
komp« 
nach | 


folgt 1 


(11) 


M. 
und 
stal 
mal 
nun 
(9) 
dv dw Dab 
(10a 
. mit 
Bolt 
Gasr 
fortp 
mise] 
wirkt 
berec 
= 
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Hierin ist p = p, + p, gleich der Summe der Partialdrucke, und c, 

und c, sind die Molkonzentrationen n,/n bzw. n,/n. Die Kon- 

stante ¥,, ist maBgebend fiir die sogenannte Thermodiffusion. Wie 

man aus (7) erkennt, ist x,, = x,, und »,, =— »,. Fiir die Span- 
nungskomponenten p,, usw. ergeben sich folgende Werte: 

ou 

2% 


| usw.; 


Dabei ist 7 der Koeffizient der inneren Reibung der Mischung. 
Für die Komponenten des Wärmestromes ergibt sich: 

oT 
Ox 
mit 2 analogen Gleichungen für q, und q,. Die Größe k ist die 
Boltzmannsche Konstante, A das thermische Leitvermögen der 
Gasmischung in Abwesenheit von Diffusion, d.h. wenn „=uw,=u. 

Damit sind sämtliche Gleichungen zur Behandlung des Schall- 
fortpflanzungsproblems gegeben. 

Enskog hat bekanntlich die Konstanten 7, /, x,,, ?,, der Gas- 
mischung unter speziellen Voraussetzungen über die Wechsel- 
wirkungskräfte zwischen den Atomen beim Stoß näherungsweise 
berechnet. 


(10a) q,= k Tin, (u, — u) +n, (u,— u] + T(u, —u,) — A 


$2. Fortpflanzung von Dichtewellen in der Gasmischung 


Wir betrachten nur den Grenzfall unendlich kleiner Amplituden, 
= (y— 1)0. Hierin ist 7 das 
Verhältnis c,/¢,, das im vorliegenden Falle einer Mischung ein- 
atomiger Gase stets °/, ist (in nicht zu tiefen Temperaturen, da 
sonst quantenmechanische Korrekturen notwendig sind). Das mittlere 
Molekulargewicht der Mischung ist M = c, M, + ¢,M,. Auch diese 
GréBe wird in der Dichtewelle wegen der auftretenden Diffusions- 
vorgänge örtlich und zeitlich schwanken. Wir setzen daher: 
M=M,(1+ u). Die Größen o, ©, u sind dimensionslos und klein 
gegen 1. 

Wir setzen in die Gl. (5a,b,c) die Werte (9) für die Spannungs- 
komponenten ein, und differenzieren die 1. Gleichung nach z, die 2. 
nach y und die 3. nach z und addieren die 3 Gleichungen. Daraus 
folgt unter Benutzung der Kontinuitätsgleichung (4): 


2 


und setzen: o= g,(1 + o) und 


(9) 
| =— + Ge] usw. 
7 
ie 
en 
(2) 
er 
de 
en 
de 
2 
die 
ngt 
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Nun ist p = se - 0, also: (v= ah , b also die Newtonsche Schall. 

geschwindigkeit) 

(12) Au). 

Hierbei sind die Produkte von Differentialquotienten von o, © und p 

vernachlässigt. Somit wird aus (11): 


(ila) — do + 40-0 


= 9 


Darin ist a = vY y die gewöhnliche Schallgeschwindigkeit, und 
Nun zur Umformung der Energiegleichung (6). Es ist: 


Q= e kTın, +n,)=oe,-T, wo c, die gewöhnliche spezifische 
Wärme bei konstantem Volumen pro Masseneinheit ist. Es wird also: 


Eu 0, T= = T, wo C, die Molwärme ist. Weiter folgt aus 
(10a), (10b), (10e): 


+ + + (¢, 0, + 0) 


oy Oz 2 
2 


OW, 
Ox 


0 ou Ov 


Oz 


Nun ist aber c, u, + ¢,u, nichts anderes als die x-Komponente u, 

der mittleren Geschwindigkeit. Nach (8) ist die Divergenz dieser 
on 

ist die Divergenz der Schwerpunktsgeschwindigkeit nach (4) gleich 


Ox 


4 dw, — v2) + | _1AT. 


dy 


Geschwindigkeit aber in unserer Näherung — 2 . Andererseits 


, somit ist die erste eckige Klammer der rechten Seite: 


n Olgm 
Moleküls ist. Damit folgt für die Divergenz des Wärmestromes: 


04; 04, 04: 
Ox + oy + Oz 


5 0 


, wo m=c¢,m,+c¢,m, die mittlere Masse eines 


Da p, = ¢,p und p, =¢,p ist, folgt für die eckige Klammer: 


Nun t 
Man 
Divisi 


(16) 


Die 3 € 
Um zu 
wir deı 
(17) 
Mit die 
(11b) 


Annak 


M.. 
Nun 
Aus 
(13) 
Also 
(14) 
: 5 | 
| 
Qı 
Für d 
Unter 
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Nun ist: 
+%,=1 und «M +«,M,=M. 


Aus diesen beiden Gleichungen folgt: 


M-M, * 


Also: 


M,— M, M 


Ap, — = % 


Unter des für 4 p ergibt sich schließlich: 


D u M, — M, 
at + at 4,40 Cı Cr 


[de 149 Au) — 


| 04 
(15) 
M, 
| 
Hierin ist gesetzt: 
A+ 12 Io » , 
= = ; Ip = 2%2 Ty. 
+ 12 = %29ı2 Ly 
Nun benötigt man noch eine 3. Beziehung zwischen den 3 Größen a, 9, u. 
Man erhält sie, wenn wir (2b) von (2a) subtrahieren nach vorheriger 
Division durch o, bzw. o;: 


6 1 


ax | 


Für die linke Seite der Gleichung kann auch geschrieben werden: 


ot @. Ot 


#@e, 


Unter TORE von (13), (7a,b,c), (8), (14) folgt: 


(16) 


| — Da du + mac - My —1)40=0. 


Die 3 Gleichungen (11a), (15) und (16) bestimmen die 3 Größen a, ©, p. 
Um zu einem Ausdruck’ für die Schallabsorption zu kommen, machen 
wir den Ansatz: 

(17) 2) = Beiwt-az), “= C 


Mit diesem Ansatz gehen wir in (11a), (15) und (16) ein, und es folgt: 


Annalen der Physik. 5. Folge. 39. 15 


h 

1 

h 

e: 
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A vj2 C2 | 
- [77 M, M 
(15a) + Bi “(y — 1) 
M, 
M, — M, M, 
— M20) C2 + |= 9: 


| 

M,. — M, \2 

M, — M, \? 

+0 |‘? + + = 0. 


Wenn diese 3 homogenen linearen Gleichungen für die 3 Ampli- 
tuden A, B, C gleichzeitig bestehen sollen, so muß die Determinante 
dieses Gleichungssystemes verschwinden. Diese Gleichung liefert 
den gesuchten Zusammenhang zwischen der Frequenz ® und der 
Wellenzahl x = 2r/Wellenlänge. Es wäre nun inkonsequent, die 
Lösung der Determinantengleichung genauer als bis zu Glieden 
1. Ordnung in »/, A, Dia, 912 zu bestimmen, da die zugrunde liegenden 
Bewegungsgleichungen (7) und (9) auch nicht genauer sind. Die 
0. Näherung erhält man, wenn sämtliche Relaxationsphänomene ver- 
schwinden, d.h. wenn 7/ = 4, = Dy = vj2 = 0. Unter diesen Voraus- 
setzungen folgt aus (16a) C = 0 und aus (11b) und (15a) ergibt sich 
4 = 1 und =a*, Auch bei Anwesenheit von Relaxations- 
vorgängen werden diese Beziehungen noch näherungsweise gültig 
sein. Wesentlich ist, daß C dann proportional zu Gliedern 1. Ord- 
nung in 7’, A, Dia, via wird. Damit folgt in genügender Naherung 
aus (16a): 


io M, — M, \? 


M, 
M,— M, M, - M, 
+ By 1 ( M, + C2 ‘|. 


Diesen Wert setzen wir in die beiden übrigen Gleichungen ei 
und erhalten . . 


| 


M, 
+ M2 C2 


M 


| 
Die 
chu 
| | (18) 
Sch 
so f 
(19) 
Koef 
i’ = 
(19a) 
Die z 
der 
sieht, 
auch 
abhär 
sehen 
| 
Diffus 
der ec 
Molko 
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uM, M,— M, 


io 
4 12% + ¢ ( 


13 
| 
i rs M. 
(15a’) | +B. + + if 
| 
| 


2 
+ Dy, - =) | =0. 


M, 
Die Determinante dieses Systems von 2 linearen homogenen Glei- 
chungen setzen wir nun gleich Null, und es folgt: 


-1 


io 


4, M, \2 
sa [guts 


'M, - 

M,- M,) 
M, 


(18) 
+ 2(7 — 1) 


Schreiben wir den Ausdruck fiir die Schallbewegung in der Form: 


2ar 
ao = Ae~m i ), 


so folgt: 


= 4, 

19) | 2a 

( ‘ 
M, — M, \2 ; M, — M, 

| + 7 ( + 2(7 — 1) 


0 


§ 3. Diskussion der Ergebnisse 


Setzt man in (19) den Diffusionskoeftizienten D,, und den 
Koeftizienten der thermischen Diffusion »j2 gleich Null, so folgt 
i, =% und es ergibt sich die Kirchhoffsche Formel: 


2 
Die zu Dy. und vig proportionalen Glieder in (19) geben den Einfluß 
der Diffusionserscheinungen auf die Schallabsorption. Wie man 
sieht, wird durch das Hinzutreten dieser Glieder die Frequenz- und 
auch die Druckabhängigkeit der Absorption nicht geändert. Bei der 
Frequenzabhängigkeit sieht man das sofort ein. Bei der Druck- 
abhängigkeit ist zu beachten, daß D,,, und wie wir später noch 
sehen werden, auch rj. dieselbe Druckabhängigkeit besitzen wie 1 
und 2’, nämlich proportional sind zu p=!. Das von der gewöhnlichen 
Diffusion herrührende Glied: 
2 

) 
der eckigen Klammer von (19) ist proportional zum Produkt aus den 
Molkonzentrationen c, und c, der Mischung und außerdem pro- 

15* 


| | 
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portional zum Quadrat der Massendifferenz. Dieses Glied ver- 
schwindet also 1. wenn es sich um ein reines Gas handelt, wie es 
sein muß (c, oder c, = 0), und 2. wenn die beiden Atomarten die- 
selbe Masse haben. Das Glied wird klein sein, wenn die Unter- 
schiede in den Massen der Atome gering sind. Da D,, von der- 
selben Größenordnung ist wie 7 und 2’, so wird die gewöhnliche 
Diffusion nur dann einen wesentlichen Beitrag zur Schallabsorption 
der Mischung ergeben, wenn die Massen der beiden Atomsorten sehr 
voneinander verschieden sind, und wenn außerdem keine der beiden 
Molkonzentrationen c, oder c, klein ist. Da der Diffusionskoeffizient D,, 
näherungsweise vom Mischungsverhältnis unabhängig ist, so erreicht 
obiges Glied seinen größten Wert wenn c, und c, ungefähr gleich 
und gleich 0,5 sind. Die restlichen in (19) auftretenden zusätzlichen 
Absorptionsglieder hängen mit der Thermodiffusion zusammen. Um 
das Verhalten dieser Glieder diskutieren zu können, müssen wir 
einen Ausdruck für »,, haben. Näherungswerte für diese Größe 
hat Enskog?) im Falle elastischer, starrer Kugeln gegeben. Sie 
lauten [unter Verwendung von (8)]: 


, nf,—n 
Vg = T = nate =)- 


(20) 
ah — af 


f 20 (m, + m,)* 
+ 16m, m, + 13mg") 16m, m, + 30m,*) 


2 


82, (m; + my) 2(m, + m,)* 
2 (30m,? + 16m,m, + 13m,*)(13m,? + 16m, m, + 30:m,?) 


| 

Vm, (m, +m,) 2 (m, + m,)* 

8; m, + my 8 (m, + my)? 

8°, 2m, + 16m, m, + 13m, ’ 

8, 8 (m, + my)? 

8: 2m, 13m,? + 16m, m; + 30m,? 

(30m,? + 16m, m, + 13m,%) (13m,? + 16m, m, + 30m) 


Hierin sind s, und s, die Durchmesser der elastischen Kugeln, die 


1) D. Enskog, Diss. S. 100ff. 
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die beiden Atomarten repräsentieren, und s,, ist das arithmetische 
Mittel von s, und s,. 

Das einzige in (19) auftretende zu vj. proportionale Glied wird 

unter Benutzung von (20): 

2.7 — Imp = 2(y — Di 
ahı = 
(Cy + + — Cy 

Die Ausdrücke für f, und f, verschwinden, wenn M, = M, und s, =s, 
ist, so daß das Glied sehr klein wird, wenn die beiden Atomarten nahezu 
gleich schwer und gleich groBsind. Es verschwindetaberauch schon, wenn 
s, +s, aber M, = M, wegen des Faktors M, — M,. Im übrigen ist das 
Glied proportional zu c,¢, und zu y—1. Daßes proportional zu 7 — I 
ist, erklärt sich daraus, daß Thermodiffusion nur bei Vorhandensein 
von Temperaturdifferenzen auftreten kann. Diese verschwinden aber, 
wenn y= 1. Wichtig ist besonders das Vorzeichen dieses Gliedes. 
Wie man sieht, ist im Spezialfall s, = s, = s,, (Fall von Isotopen): 
j>0, fo< 0, wenn das Gas 1 das schwerere ist. Unter Voraus- 
setzung elastischer Kugeln und gleich großer Atome ist also via 
positiv, wenn M, > M, ist. Damit wird unser betrachtetes Glied 
negativ. Die gewöhnliche Diffusion und die thermische Diffusion 
wirken in diesem Fall einander entgegen. Im übrigen ist dieses 
Glied von ähnlicher Bauart wie das von der gewöhnlichen Diffusion 
herrührende Glied. Im Falle der Potenzgesetze hat Enskog’) 
gezeigt, daB für gleich große Atome und M, > M, die Größe »,, 
positiv ist, wenn der Exponent des AbstoBungsgesetzes größer als 5 
ist (darin ist auch der Fall elastischer Kugeln enthalten), und negativ 
wird, wenn dieser Exponent kleiner als 5 ist. Im Falle eines Ab- 
stoBungsexponenten kleiner als 5 ist dann unser betrachtetes Glied 
positiv, und gewöhnliche und thermische Diffusion wirken im 
gleichen Sinn. 

Wir wollen uns noch klar machen, warum gewöhnliche und 
thermische Diffusion im entgegengesetzten oder im gleichen Sinne 
wirken, je nachdem rj. positiv oder negativ ist. Dies erkennt man 
leicht an Hand der Diffusionsgleichung (7a). Führen wir p, = ¢,p 
und p, = c,p ein, und beachten, daß die Gradienten von c, und c, 
verhältnismäßig klein sind, so folgt: 


— QC, Op 


s 


| dz 


1) D. Enskog, Arkiv f. Mat. Astr. Fysik 16. 8. 55. 1922. 


mm, Ox 


. 
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In einer Dichtewelle hat nun grad p in 1. Näherung dieselbe Richtung 
wie grad T, solange der Vorgang rein adiabatisch verläuft. Da 
der Faktor von Op/Oz in der Klammer im Falle M, > M, stets 
negativ ist, so folgt, daß gewöhnliche Diffusion und Thermodiffusion 
entgegengesetzt wirksam sind, wenn »,, positiv ist, und im gleichen 
Sinne wirken, wenn r,, negativ ist. Da «,, stets positiv ist, folgt 
im 1. Fall vj, > 0, und im 2. Fall xj. < 0. Die gewöhnliche Diffusion 
hat zur Folge, daß die schwereren Teilchen sich in Richtung zu- 
nehmenden Druckes bewegen. 

Nun zur Diskussion des dritten in der eckigen Klammer von (19) 
auftretenden zusätzlichen Absorptionsgliedes. Es lautet: 


A 


7 Yo Cr 7 

Es ist zunächst klar, daß dieses Glied stets positiv ist. Im Falle 
elastischer Kugeln erhält man: 

Dieses Glied verschwindet, wenn y = 1, oder wenn c, oder c, =0, 
oder wenn s,=s, und M, = M, ist. Es verschwindet dagegen 
nicht, wenn s, =s, aber M, + M,, oder wenn M, = M, aber 
8, + 8,. Es ist im allgemeinen von derselben Größenordnung wie die 
beiden schon diskutierten Glieder. Daß es für y = 1 verschwindet 
hängt damit zusammen, daß Thermodiffusion in diesem Falle gar 
nicht auftreten kann. 

Die Existenz dieses Gliedes ist auf eine durch Thermodiffusion 
vergrößerte Wärmeleitfähigkeit zurückzuführen. 

Bisher wurden bei der Diskussion der einzelnen Glieder von (19) 
die Enskogschen Näherungsformeln für den Fall elastischer Kugeln 
benutzt. In einem 2. Fall lassen sich sogar exakte Werte für den 
Koeffizienten »,, angeben, und zwar in dem Grenzfalle, wo M,/M, 
und s,/s, klein gegen 1 sind. Dieser Grenzfall entspricht demjenigen, 
der in der Lorentzschen Form der Elektronentheorie der Metalle 
vorliegt. Bedeutet » den Exponenten des Wechselwirkungsgesetzes, 
so ist ): 

(21) v 


36-1) 

Damit folgt fiir die 3 durch Diffusion bedingten Glieder der 
eckigen Klammer von (19) der Reihe nach, wie sie bisher diskutiert 
wurden: 


nik. 


1) D. Enskog, Arkiv f. Mat. Astr. Fysik 16. 8.56 und Ann. d. Phys. 38. 
S. 750. 1912. 
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7+ Diy (gewöhnliche Diffusion), 


- Ne, D (Thermodiffusion), 


12 
p (bedingt durch vergrößerte Wärmeleitung 
infolge thermischer Diffusion). 
Beim 3. Glied tritt im Nenner ein c, auf. Für c, = 0 muß die 
Annäherung auch tatsächlich versagen, da die Verteilungsfunktion 
der leichten Atome nur durch Zusammenstöße mit den schweren 
Atomen hergestellt wird, was natürlich nicht mehr möglich ist, 
wenn c, =0 ist, d. h. wenn keine schweren Atome mehr vor- 
handen sind. 


In einem weiteren Spezialfall lassen sich noch strenge Ergebnisse 
erhalten, nämlich wenn sich die Atome gemäß der 5. Potenz ihrer 
Entfernung abstoßen. Dann ist bekanntlich »,, = 0, und es bleibt 
in der eckigen Klammer von (19) nur das gewöhnliche Diffusions- 
glied übrig, da in diesem Fall die Thermodiffusion verschwindet. 


Besondere Betrachtung verdienen noch die Isotopengemische. 
Da bei diesen im allgemeinen die Massen nur wenig voneinander 


| mM, 


‘verschieden sind, d. h. da <1, so folgt aus unseren 


Betrachtungen im Falle PR... Kugeln, daß der Beitrag der 

Diffusionserscheinungen zur Schallabsorption meist sehr klein sein 

dürfte. Da bei Isotopen s, =s, ist, so folgt aus den früheren 

Formeln, daß der Beitrag der Diffusion zur Schallabsorption pro- 

portional zu A wird. Ausnahmen könnten nur auftreten, 
0 


wenn es sich um Isotopengemische niedriger Ordnungszahl handelt, 


oe von der Größenordnung 1 wire. 
‘ 

Wenn man die fiir einatomige Gasgemische gefundenen Formeln 
auch auf Mischungen zweiatomiger Molekiile SUR will, z. B. 


a a | ungefähr gleich = , und das von der ge- 


wo dann | 


auf Luft, so wäre 
wöhnlichen Diffusion herrührende Glied in 1) wird dann durch den 


Faktor I 80 klein, so daß es nur noch ganz geringe 


A 50 
Beiträge zur Gesamtabsorption liefern kann. Der von Rocard 
angegebene Wert, wonach die Diffusion im Falle von Luft einen 
Beitrag liefert, der etwa gleich groß ist, wie derjenige, der von 
der Wärmeleitung herrührt, ist demnach viel zu groß. 
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Nun noch ein Vergleich zwischen unseren allgemeinen Ergebnissen 
und denen von Rocard. Der DiffusionseinfluB ist nach Rocard 
gegeben durch: 


m = AD 


12 9453 
i 


wo A= “et . Wie man sieht, stimmt die Frequenzabhängigkeit 
und Druckabhängigkeit mit der unserigen überein. Im übrigen 
kann diese Formel nicht richtig sein, da sie z.B. schon die Ab- 
hängigkeit von der Konzentration nicht richtig wiedergibt. Sie 
ergibt z. B. auch für das reine Gas (c, = 1) einen Diffusionseinfluß, 
Auch sonst gibt sie die Verhältnisse unzureichend und unzutreffend 
wieder, wie wir gesehen haben. 2 


Zusammenfassung 


Der Eintluß der gewöhnlichen Diffusion und der Thermodiffusion 
auf die Schallabsorption in einatomigen Gasgemischen wird unter 
Verwendung der Enskogschen Gleichungen untersucht. Es wird 
gezeigt, daß die Rocardschen Untersuchungen über diesen Gegen- 
stand unzureichend und nicht fehlerfrei sind. Im $ 1 werden die 
Bewegungsgleichungen für die einatomige Gasmischung auf Grund 
der kinetischen Gastheorie abgeleitet. Es ergeben sich 9 a 
Gleichungen für die 9 Unbekannten u,, v,, W,; U,, %y, W,; 9, und 9, 
und T. Im $ 2 wird dann die Ausbreitung von Dichtewellen 
unendlich kleiner Amplitude betrachtet. Eine Schallwelle in der 
Mischung äußert sich nicht nur in periodischen Schwankungen der 
Dichte und der Temperatur, sondern auch in Schwankungen des 
mittleren Molekulargewichtes derMischung. DasSystem der Bewegungs- 
gleichungen läßt sich auf ein System von 3 linearen homogenen 
Gleichungen für die Dichte, die Temperatur und das mittlere Mole- 
kulargewicht reduzieren. Die Determinantengleichung desjenigen 
Gleichungssystems, das aus dem soeben genannten durch den Ansatz 
ebener Schallwellen folgt, gibt die gesuchte Absorption. Zu den von 
Kirchhoff angegebenen Gliedern, die von der Reibung und der 
Wärmeleitung herrühren, treten 3 neue Glieder, die dieselbe Frequenz- 
und Druckabhängigkeit besitzen, wie die Kirchhoffschen Glieder. 
Das eine Glied trägt dem Einfluß der gewöhnlichen Diffusion Rechnung. 
(> — M, 


Es ist proportional zum Diffusionskoeffizienten D,,, zu M 
0 


(M = Molekulargewichte), zum Verhältnis c,/c, der spez. Wärmen 
und zum Produkt c,c, der Molkonzentrationen. Das 2. Glied trägt 
dem Einfluß der thermischen Diffusion Rechnung. Es ist proportional 
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M,-M 
mu " —1, zu ¢,¢,, zu u ' und zu einem Koeffizienten vj, der 


die thermische Diffusion bestimmt. Das Verhalten dieses Gliedes ist 
untersucht unter der Annahme elastischer, starrer Kugeln und außer- 
dem unter der Annahme, daß die eine Molekülart unendlich klein 
und unendlich leicht gegenüber der anderen ist. Wenn v7 > 0 
und M, > M, ist, dann wird das betrachtete Glied negativ’ Gewöhn- 
liche und thermische Diffusion wirken einander entgegen. Im 
Falle vj, < 0 wirken gewöhnliche Diffusion und thermische Diffusion 
im gleichen Sinn. Die gewöhnliche Diffusion bewegt die schwereren 
Teilchen in Richtung zunehmenden Druckes. Das 3. Glied ist stets 
positiv, und trägt der durch Thermodiffusion vergrößerten Wärme- 
leitung Rechnung. Auch das Benehmen dieses Gliedes ist im Spezial- 
fall elastischer Kugeln und in dem Grenztall, daß die eine Molekülart 
unendlich klein und unendlich leicht gegenüber der anderen ist, 
diskutiert. Die gefundenen Schallabsorptionseffekte fallen nur dann 
ins Gewicht, wenn die Massen oder Molekülgrößen sehr verschieden 
sind. Im Falle von Isotopengemischen ist ihr Beitrag zur Schall- 
absorption im allgemeinen nur sehr gering. Dasselbe dürfte für 
Luft gelten. Prinzipiell hat man die Möglichkeit, durch Schall- 
absorptionsmessungen an Gasmischungen den sonst der Messung 
nur schwer zugänglichen Koeffizienten der Thermodiffusion zu finden. 


Berlin, 1. Institut für theoretische Physik der Universität. 


(Eingegangen 10. Januar 1941) 
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Keimbildung in wäßriger Lösung 
Von U. Dehlinger und E. Wertz 
(Aus dem II. physikalischen Institut der Technischen Hochschule in Stuttgart) 


(Mit 4 Abbildungen) 


Die Keimbildung in Dämpfen und in kristallinen Stoffen wurde 
in den letzten Jahren vielfach studiert, dagegen sind über spontane 
Keimbildung in übersättigten wäßrigen Lösungen keine neueren Ar- 
beiten bekannt geworden. Im folgenden soll für Alaunlösungen ge- 
zeigt werden, daß durch systematische Reinigung der Lösung von 
Fremdkeimen sehr gut definierte Ergebnisse erhalten werden können, 
die einfach theoretisch auswertbar sind und die allgemeine Frage 
der Keimbildung in kondensierten Stoffen ganz neu beleuchten. 

Die Arbeiten über wäßrige Lösungen bis zum Jahre 1902 wurden 
von Wilhelm Ostwald!) zusammenfassend bearbeitet. Er kommt 
zu dem Schluß, daß in dem Temperaturkonzentrationsgebiet der 
Übersättigung eine Grenzkurve besteht, so daß bei kleineren Konzen- 
trationen keine spontane Keimbildung eintritt. Er unterscheidet also 
stabile und metastabile übersättigte Lösungen. Die Frage, wie scharf 
und für wie lange Wartezeiten diese Grenzkurve definiert ist, läßt 
er ausdrücklich offen. Auch später wurde diese Frage nicht mehr 
experimentell untersucht’). Der einzige, der überhaupt die 
Keimbildung nach längeren Wartezeiten untersucht hat, war 
de Coppet [1875°),. Er findet, daß die mittlere Lebensdauer einer 
übersättigten Glaubersalz- und Natriumbromidlösung um so kleiner 
ist, je größer die Konzentration und außerdem, je größer das Volumen 
des Gefäßes ist, wie es für einen statistischen Vorgang zu erwarten ist. 

Wir stellten uns die Aufgabe, die Existenz der Ostwaldschen 
Keimbildungsgrenze nach langen Wartezeiten sowie ihre Temperatur- 
abhängigkeit zu untersuchen. Wir verwandten Alaunlösungen wegen 
der bekannten starken „Kristallisationskraft“ dieses Salzes, welche 
die einmal gebildeten Keime sehr schnell heranwachsen und da- 
durch erkennen läßt. 


1) W. Ostwald, Allgem. Chem. II, 2; S. 752. 1920. 2. Aufl. 

2) Vgl. die neuere Zusammenfassung bei A.E.H. Tutton. Chrystallo- 
graphy and Practical Crystal Measurement. 1911. 

3) M. de Coppet, Ann. de Chim. Ser. 5; 6. S. 275—288. 1875. 
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I. Vorversuche 


Es wurde zunächst an zalılreichen Lösungen von verschiedenem 
Übersättigungsgrad festgestellt, welche Wartezeiten bei Zimmer- 
temperatur für eine spontane Keimbildung notwendigsind. Die Lösungen 
wurden aus reinstem, kristallisiertem Kalialaun (KAISO,), - 12H,0), 
der von zwei verschiedenen Firmen bezogen war, bei 80°C her- 
gestellt und auf 21° © abgekühlt. Für die Beobachtungen stand 
ein Kellerraum zur Verfügung, der lange Zeit auf 21 + 1° gehalten 
werden konnte. Die Lösungen wurden in Gläsern von je 25 ccm 
Inhalt aufbewahrt, die mit Gummidichtungen abgeschlossen waren. 
Die Bildung von Keimen konnte durch das Glas hindurch leicht 
festgestellt werden. 

Anfangs erhielt man für gleiche Proben stark schwankende 
Wartezeiten; auch verhielt sich das von den beiden Bezugsquellen 
stammende Material nicht gleichartig. Darauf wurden die heißen 
Lösungen filtriert, was aber die Ergebnisse kaum änderte. Allmäh- 
lich fand man aber, daß die Wartezeiten wesentlich länger und 
besser reproduzierbar werden, wenn man die Lösungen bei der Her- 
stellung längere Zeit auf 80° hält oder wenn man die Kristalle 
mehrmals umkristallisiert. Z.B. zeigte eine Lösung, die bei 21° 
eine Übersättigung von 170°/, besaß!) und die bei der Herstellung 
20 Min. auf 80° gehalten war, eine Lebensdauer von 87 Min., da- 
gegen hatte dieselbe Lösung, 4 Std. bei 80° gehalten, eine Lebens- 
dauer von 60 Std. Wird die durch kurzzeitiges Erhitzen hergestellte 
Lösung nach dem Auskristallisieren wieder nur 10 Min. auf 80° ge- 
bracht und dadurch wieder gelöst, so verlängert sich die Lebens- 
dauer von 87 Min. auf 7—20 Std. 

Bei dem langdauernden Erhitzen auf 80°, noch mehr bei höheren 
Temperaturen, entsteht eine leichte Trübung, die auf hydrolytische 
Spaltung des Alauns und kolloide Ausfällung der Spaltungsprodukte 
zurückzuführen ist. Es erhob sich die Frage, ob die Verlängerung 
der Wartezeit beim Erhitzen nicht darauf zurückzuführen ist, daß 
diese Kolloide oder andere Veränderungen der Lösung die spontane 
Keimbildung hemmen. Dies wurde dadurch widerlegt, daß Lösungen, 
die durch Erhitzen eine lange Wartezeit erhalten hatten, mit solchen 
von kurzer Wartezeit gemischt wurden. Es ergab sich stets eine 
kurze Wartezeit für die Mischung. Es befindet sich also nicht etwa 
in den lang erhitzten Lösungen ein die spotane Keimbildung 
hemmendes Agens, sondern umgekehrt sind in den ursprünglichen 
Lösungen keimbildende Verunreinigungen, welche durch mehrmaliges 


1) Nach den Angaben über Löslichkeit in Landolt-Börnstein Bd. I. 8. 632. 
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Umkristallisieren sowie durch langes Erhitzen der Lösung zerstört 
oder entfernt werden können. 

In kleinen Mengen zugefügte kolloide Stoffe, z. B. Gummi, 
AIOH),, geben keine merkbare Veränderung der Wartezeiten. Da. 
gegen führen kristallisierte Stoffe, z. B. NaClO, schnell zur Kristallisation, 
Im übrigen haben wir uns mit der Frage nach der Natur der in 
den ursprünglichen Lösungen vorhandenen, vermutlich von der Her. 
stellung des Alauns herrührenden Verunreinigungen nicht weiter 
beschäftigt, sondern die in den umkristallisierten und genügend 
lang erhitzten Lösungen auftretende spontane Keimbildung näher 
untersucht. 


Auf Grund dieser Vorversuche wurden für die Hauptversuche 
Lösungen aus Kristallen benützt, die durch 3maliges Umkristallisieren 
aus heißer Lösung gewonnen waren. Durch besondere Versuche 
wurde festgestellt, daß ein weiteres Umkristallisieren die Wartezeit 
nicht mehr merkbar verlängert. Im allgemeinen wurden die Lösungen 
jeweils heiß filtriert. Nur bei den höchsten Konzeutrationen war 
dies nicht mehr möglich, da sich hier infolge Verdampfung während 
des Filtrierens schon spontan Keime bildeten; jedoch haben die 
Vorversuche gezeigt, daß das Filtrieren die Wartezeiten nicht merk- 
bar ändert. 


Mit den so hergestellten Lösungen erhielt man die im folgenden 
Abschnitt dargestellten reproduzierbaren Ergebnisse. Nur 3 von 
insgesamt 154 Proben verhielten sich auch nach der Reinigung ab- 
weichend. Sie enthielten Lösungen mit 64, 72 und 180°/, Über- 
sättigung und bildeten nach 134 bzw. 128 und 237 Tagen Kristalle, 
während alle andern Proben gleicher oder ähnlicher Konzentration 
bis jetzt, nach mehr als einem Jahr, überhaupt nicht kristallisierten. 
Es ist anzunehmen, daß diese ganz aus der sonstigen Gesetzmäßig- 
keit herausfallenden Proben Fremdkeime enthielten, die zufällig auch 
bei dem wiederholten Umkristallisieren nicht entfernt wurden. 


Il. Feststellung der metastabilen Grenzkonzentration bei 21°C 


Es wurden Lösungen mit einer Übersättigung von 4—620), 
hergestellt und so wie oben angegeben, von keimbildenden Ver 
unreinigungen befreit. Bei 21° erhielt man mit ihnen die in Tab. | 
angegebenen Wartezeiten bis zur Keimbildung. Lösungen mit einer 
Konzentration von weniger als 3,83 Mol.-Prozent (520°/, Übersättigung 
bildeten mit den am Schluß von Abschnitt I erwähnten, deutlich 
herausfallenden Ausnahmen überhaupt keine Keime. Ein große 
Teil dieser Lösungen steht jetzt mehr als 1 Jahr unter Beobachtung 


Abb. 
bildur 
21°C. 

mit 


bei e 
warte 
Kurve 
die ] 
Losur 
daher 
durch 
zuges« 
sponte 
Grenz 
iibersi 
der er 


Wa 
geh 
der 
Wie 
Pun 
| sich 
Asy 
um 
send 
| 


U. Dehlinger u. E. Wertz. Keimbildung in wäßriger Lösung 229 


In Abb. 1 sind die Mittelwerte der 
Wartezeiten von Tab. 1 mit den zu- 
gehörigen Schwankungen als Funktion 
der molaren Konzentration aufgetragen. 
Wie man sieht, läßt sich durch die 
Punkte eine glatte Kurve durchlegen, die 
sich einer bei 3,83 Mol.-Prozent liegenden 
Asymptote nähert. Die Schwankungen 
um die Mittelwerte nehmen mit wach- 
sender Steilheit der Kurve zu, wie es 
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Abb.1. Mittlere Wartezeiten bis zur Keim- 

bildung in iibersittigten Alaunlösungen bei 

21°C. Ordinate: Mittelwert der Wartezeit 
mit Angabe des lfach mittleren Fehlers 


bei einem statistischen Vorgang zu er- 
warten ist. Die Existenz der glatten 
Kurve für die Mittelwerte beweist, daß 
die Keimbildung bei den gereinigten 
Lösungen gesetzmäßig verläuft; sie kann 
daher nicht zufälligen Einwirkungen 
durch Verunreinigungen und ähnliches 
zugeschrieben werden, sondern muB 
spontan sein. Die Asymptote gibt die 
Grenze des metastabilen Gebiets der 
übersättigten Lösungen. In Anbetracht 
der erwähnten Beobachtungszeit von mehr 
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als einem Jahr kann gesagt werden, daß diese Grenze fiir die Keim. 
bildung nach den bisherigen Erfahrungen ebenso scharf definiert ist 
wie die Löslichkeitslinie selbst. Es ist daher zu vermuten, daß 
auch sie durch eine bestimmte thermodynamische Bedingung gekenp. 
zeichnet ist, auf die weiter unten eingegangen wird. 


Die Auskristallisation nach der spontanen Keimbildung erfolgte 
in Anbetracht der hohen Konzentrationen sehr schnell, so dab 
manchmal sogar das Gefäß gesprengt wurde. Sie ergab vielkörnige, 
knollige Kristallaggregate, wie bei Anwesenheit zahlreicher Keime zu 
erwarten ist. Dagegen ergaben die oben beschriebenen heraus- 
fallenden Proben, wie auch die sonstigen nicht durch spontane 
Keimbildung kristallisierten Lösungen einzelne klare Kristalle, wie 
sie für wenige Keime kennzeichnend sind. 


Die Freundlichkeit der Herren Gerhard Schmid und 
U. Jetter hat es uns ermöglicht, einige nicht kristallisierende, über- 
sättigte Lösungen in einem intensiven Ultraschallfeld zu untersuchen. 
Auch hier — wie übrigens auch beim gewöhnlichen Schütteln — 
ergab sich keine spontane Kristallisation. Dagegen beschleunigte 
das Ultraschallfeld die Kristallisation an einem von außen herein- 
gebrachten Keim stark, und zwar besonders dadurch, daß es die 
sich bildenden Kristalle in kleine Stücke aufteilte und in der Lösung 
zerstreute. Diese Versuche zeigen die hohe Stabilität der über- 
sättigten Lösungen unterhalb der metastabilen Grenze besonders 
deutlich. 


III. Die Temperaturabhängigkeit der metastabilen Grenze 


Nachdem die Existenz und Lage der Grenze für die spontane 
Keimbildung bei 21° festgestellt war, wurden zur Messung ihrer 
Temperaturabhängigkeit die Wartezeiten bei 0 und 40°C bestimmt. 
Die ebenso wie früher vorbehandelten Proben wurden in einen Eis 
thermostaten und in ein durch einen Regler auf 40 + 0,5° C ge 
haltenes Wasserbad gebracht. Die erhaltenen Wartezeiten sind in 
Tab.2, ihre Mittelwerte in Abb.2 eingetragen. Da weniger Proben 
untersucht werden konnten als bei Zimmertemperatur, ist die Ge 
nauigkeit der Grenzbestimmung nicht ganz so groß wie dort. Jedoch 
kann folgendes gesagt werden: Während für 21° die Grenze sicher 
zwischen 3,7 und 3,83 Mol-Prozenten KAI(SO,), liegt, liegt sie bei 
40°C sicher unter 4,0 Mol-Prozenten; bei 0° C liegt sie wenig unter 
3,5 Mol-Prozenten. Aus Abb.3 ist zu ersehen, wie sich danach die 
Temperaturabhängigkeit der metastabilen Grenze im Vergleich zur 
Löslichkeitsgrenze gestaltet. 
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Tabelle 2 


Wartezeiten bis zur Keimbildung in Stunden 


E _Mittel- 
Mol-"/, Bei 0°C werte | Bei 40°C 
KAI(SO,), in Std. | 


6 6 
22 46 
69) 94, 92, 
192! 228) 216216 204 
> 250 > 250) > 250 


§ § § § R 


Mol- % Mol X HAl (Sl 
Abb. 2. Wartezeiten bis zur Abb.3. Experimenteller Verlauf der 
Keimbildung in übersättigten Löslichkeitskurve und der metastabi- 
Alaunlösungen bei 0 und len Grenzkurve. Gestrichelt der 
40°C theoretisch bei höheren Temperaturen 
zu erwartende Verlauf der letzteren 

Kurve 


IV. Anwendung der Keimbildungstheorie nach Volmer 


Für die Keimbildung in übersättigtem Dampf liegt auf Grund 
der Arbeiten von Gibbs, Haber, Volmer, R. Becker und Döring 
u.a. eine ausgearbeitete thermodynamisch -statistische Theorie vor, 
über die in letzter Zeit A. Volmer!) ausführlich berichtet hat. Er 
betont dabei mehrfach), daß diese Theorie ohne irgendein hypotheti- 
sches Moment ableitbar und daher auf alle Fälle von Keimbildung 
in ihren Grundzügen anwendbar sei. Dies ist nicht richtig‘); so 


1) A. Volmer, Kinetik der Phasenbildung, Dresden und Leipzig 1939. 

2) Vgl. S. 20 und $. 214 des angeführten Buchs. 

3) Trotzdem wurden schon von verschiedenen Seiten die Ansätze Vol- 
mers ohne weitere Diskussion als Grundlage für theoretische Überlegungen 
im Gebiet des kristallinen Zustands benützt. Vgl. G. Masing, Ztschr. f. 
Metallkde. 32. S. 19. 1940. 
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haben die experimentellen Untersuchungen über die Kaltaushärtung 
in metallischen Mischkristallen auf eine ganz andere, in dem Buch 
von Volmer überhaupt nicht erwähnte Art von Keimbildung ge. 
führt, deren thermodynamische Grundlagen zuerst von U.Dehlinger!), 
später von R. Becker?) dargestellt wurden. 

Zwischen den beiden genannten Keimbildungsmechanismen be- 
steht folgender grundsätzliche thermodynamische Unterschied: Bei 
der Keimbildung in Dämpfen entsteht ein Keim durch Zusammen- 
stéBe von Molekülen in vollständig kompakter Form. Sein thermo- 
dynamisches Potential ist daher gleich dem der vollständigen neuen 
Phase, vermehrt um den Beitrag der Oberflächenenergie. Die thermo- 
dynamische Stabilität des Keims wird also nur bestimmt durch seine 
Größe und den Übersättigungsgrad. 

Demgegenüber entsteht ein Keim nach Dehlinger dadurch, 
daß oberhalb einer bestimmten Übersättigungsgrenze, wie man exakt 
thermodynamisch nachweisen kann, der Diffusionskoeffizient negativ 
wird. Daher werden jetzt kleine Konzentrationsschwankungen, die 
statistisch immer vorhanden sind, allmählich anwachsen. Der Keim 
bleibt also stetig mit der übersättigten Phase verbunden°); er ist 
nicht durch einen Konzentrationssprung abgegrenzt, wie der Keim 
Volmers, sondern die Konzentration eins oder nahezu eins im 
Innern des Keims geht allmählich über in die des Mischkristalls. 
Infolgedessen wird die thermodynamische Stabilität des Keims auch 
durch das negative Entropieglied erhöht, welches dieser von Null 
und eins verschiedenen Konzentration zukommt. Daher ist der Keim 
oberhalb der erwähnten Übersättigungsgrenze bei jeder Größe stabil, 
er braucht also nicht wie der Keim Volmers durch Zusammen- 
stéBe, von Molekülen innerhalb einer statistischen Schwankungs- 
periode gebildet zu werden, sondern kann durch einzelne Atom- 
sprünge allmählich aufgebaut werden. 

Kinetisch entstehen die Keime Volmers dadurch, daß einzelne 
Moleküle nach Größe und Richtung außergewöhnliche Geschwindig- 
keiten besitzen. Demgegenüber brauchen bei der Keimbildung nach 
Dehlinger die einzelnen Bewegungsschritte bei der Diffusion nur 


1) U. Dehlinger, Ztschr. f. Physik 102. S. 633. 1936; Ztschr. phys. 
chem. Unterr. 50. 134. 1937; Ztschr. f. Metallkde. 29. S. 401. 1937. 

2) R.Becker, Ztschr. f. Metallkde. 29. S.345. 1937; vgl. auch G. Borelius, 
Ann. d. Phys. [5] 28. S. 507. 1937, wo aber der Zusammenhang mit der Dif- 
fusion nicht erkannt ist. 

3) Diese stetige Verbundenheit zeigen bei der Kaltaushärtung vor allem 
die Röntgenaufnahmen, vgl. U. Dehlinger, Chemische Physik der Metalle 
und Legierungen, Leipzig 1939. 
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iber Strecken von der Größe eines Atomabstands zu führen; größere 
Schwankungen der Geschwindigkeit, die Atomsprünge über mehrere 
Atomabstände (Potentialschwellen) ergeben würden, würden den Keim 
in seinen Anfangsstadien wieder zerstören. Die Keimbildung nach 
Dehlinger wird also dann vor der von Volmer bevorzugt sein, 
wenn die statistischen Schwankungen der Geschwindigkeit klein sind 
gegenüber denjenigen der Lage, so wie es in kondensierten Stoffen 
allgemein der Fall ist. Es ist daher zu vermuten, daß in konden- 
sierten Stoffen überhaupt keine Keimbildung nach Volmer vor- 
kommt ). 

Wir untersuchen nun die Frage, ob unsere bei der Alaunlösung 
gewonnenen experimentellen Ergebnisse mit der Theorie Volmers 
in Übereinstimmung zu bringen sind. Dagegen spricht zunächst 
schon die Schärfe der von uns gefundenen metastabilen Grenze, da 
die Keimbildungstheorie Volmers, wie er ausdrücklich betont, keine 
ganz scharfe Grenze erwarten läßt. 

Eine sichere Entscheidung ist möglich, wenn wir die Tem- 
peraturabhängigkeit der metastabilen Keimbildungsgrenze mit der 
der Löslichkeitsgrenze vergleichen. Die letztere ist thermodynamisch 
durch folgende Gleichung zwischen dem thermodynamischen Potential 
der Lösung und des Kristalls bestimmt: 


Für das thermodynamische Potential des Kristalls schreiben wir: 


86, _ 
92° 


Für die Lösung schreiben wir in bekannter Weise 


= 9, +RTInc, 


wo c die Konzentration in Molen ist. 


1) Von A. Volmer (a.a.O.) werden Messungen der Temperaturabhängig- 
keit der Ausscheidungsgeschwindigkeit bei Au-Pt nach der Widerstands- 
methode von C.H. Johansson und C. Hagsten, Ann. d. Phys. [5] 28. S. 520. 
1937, mit seiner Keimbildungstheorie ausgewertet. Diese Messungen enthalten 
aber keinen Beweis dafür, daß der Volmersche Keimbildungsmechanismus 
in diesem Fall überhaupt vorliegt; vermutlich erfassen sie die Geschwindig- 
keit des Keimwachstums mehr als die der Keimbildung. Aufklärung über 
die Ausscheidungsvorgänge in Legierungen kann man nur durch Kombination 
mehrerer Beobachtungsmethoden erhalten. Die ausgedehnte Literatur auf 
sew Gebiet wird aber von Johansson und Hagsten iiberhaupt nicht 
erwähnt. 


Annalen der Physik. 5. Folge. 39. 16 


8G, 

On On 

| | 
> 


234 Annalen der Physik. 5. Folge. Band 39. 1941 
Damit wird die Gleichung der Löslichkeitsgrenze 
RTlnc =9, —9, 


Wenn wir, wie üblich, die Volumänderung beim Kristallisieren yer. 
nachlässigen, also thermodynamisches Potential und freie Energie 
vertauschen, bedeutet die Größe q die partielle molare Lösungs- 
wärme. Bei strenger Gültigkeit der Neumann-Koppschen Regel 
von der Additivität der spezifischen Wärmen wäre sie von der Tem- 
peratur unabhängig. Wir können also damit rechnen, daß q wenig 
temperaturabhängig ist. Für ideale Lösungen wäre g, und damit 
auch q von der Konzentration c unabhängig. Bei nicht idealen 
Lösungen wird es somit wenig von c abhängen. In dem Ausdruck 
für die Neigung der Löslichkeitskurve 
RT, 0a 
c dc 


Rine- 


sind daher 3 und 2% kleine Größen. 
c 07 
Das Gleichgewicht zwischen der Lösung und einem Keim, 

dessen thermodynamisches Potential sich von dem des kompakten 
Kristalls um die (auf ein Mol umgerechnete) Oberflächenenergie 0 
unterscheidet, wird dann bestimmt durch die Gleichung 

dG, 06, 

Also 


RTinc=—q+0. 


Für = ergibt sich genau derselbe Ausdruck mit den gleichen 


Größen wie oben, da O von der Temperatur und Konzentration 
nicht abhängt. 

In der Theorie Volmers ist die metastabile Grenzkurve auf- 
zufassen als die Gleichgewichtskurve derjenigen Keime, die gerade 
noch mit merkbarer Geschwindigkeit wachsen können. Daher gibt 
die obige Gleichung den nach Volmer zu erwartenden Verlauf der 
metastabilen Grenzkurve wieder. 

Wie man aus Abb. 3 leicht sieht, ist der dort dargestellte ex- 
perimentelle Verlauf der Löslichkeits- und metastabilen Grenzkurve 
mit dieser Theorie nicht zu vereinbaren. Insbesondere ist die Theorie 
nicht imstande, die Tatsache wiederzugeben, daß unterhalb c = 4°], 
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die letztere Kurve wesentlich steiler als die erstere verläuft. Wären 
aq 
de und aT gleich null, so wire = 


-. Dann müßte die metastabile Grenzkurve bei gleichem c stets 


die beiden kleinen Größen 


facher verlaufen als die Löslichkeitskurve. Wenn sie in Wirklichkeit 
wesentlich steiler verläuft, so kann dies nicht von dem Einfluß der 


kleinen Größen sf und - 


oT 
Theorie aufgegeben werden. Es ist somit festzustellen, daB sowohl 
die experimentell gefundene Schärfe der metastabilen Grenzkurve 
wie auch ihre Temperaturabhängigkeit durch die Keimbildungstheorie 
nach Volmer auch nicht annähernd wiedergegeben werden können. 


herrühren, sondern es muß die ganze 


V. Theorie der Keimbildung durch negative Diffusion 


Um eine brauchbare Theorie für die spontane Keimbildung in 
den wäßrigen Alaunlösungen zu erhalten, gehen wir von folgenden 
experimentellen Tatsachen aus: 1. Die Grenzkonzentration der 
spontanen Keimbildung ist thermodynamisch scharf bestimmt und 
von der Temperatur wenig abhängig. 2. Sieliegt wenig unterhalb des 
Molenbruchs C = 4°/,, das ist eine Konzentration, bei welcher gerade 
jedes der vier Ionen, in die das Alaun KAI(SO,), in der Lösung zer- 
spalten ist, von sechs Wassermolekülen umgeben ist. Da auch im 
kristallisierten Alaun jedes K- und Al-Ion von 6H,O umgeben ist’), 
können wir annehmen, daß diese 6 Wassermoleküle eine ge- 
schlossene Hülle um jedes Ion in der Lösung bilden. 3. Die Keim- 
bildungsgrenze nähert sich mit wachsender Temperatur der so be- 
stimmten Konzentration immer mehr. 


Wie das Experiment ergibt, ist also eine Keimbildung erst bei 
Konzentrationen möglich, die so hoch sind, daß nicht mehr jedes 
gelöste Ion seine eigene Wasserhülle besitzt, mit anderen Worten 
erst dann, wenn einzelne lonen der Lösung nur noch durch eine 
Schicht von Wassermolekülen getrennt sind. Der Keimbildungsvor- 
gang ist somit in erster Linie von der örtlichen Verteilung der 
gelösten Ionen, d. h. von der Zahl der Potentialschwellen abhängig, 
welche die Ionen bis zu ihrer Vereinigung überwinden müssen, 
während nach Ausweis der Temperaturabhängigkeit die Besonder- 
heiten der Verteilung der kinetischen Energie keine wesentliche 
Rolle spielen. Nun wird ganz allgemein die Konzentrationsveränderung, 


1) Vgl. Strukturbericht der Ztschr. f. Kristallogr. 3. S. 208, 1937. 
16* 
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die durch das fortgesetzte Uberschreiten einzelner Potentialschwellen 
entsteht, von der Diffusionsgleichung 


dc 
D 


bestimmt. Dabei ist der Diffusionskoeffizient D im allgemeinen von 
der Konzentration c abhängig, nur bei idealen Lösungen ist er von c 
unabhängig. Wir haben also zu fragen, unter welchen Umständen 
diese Differentialgleichung zu Keimbildung, d.h. zu einem allmählichen 
Ansteigen der Konzentration an einzelnen Stellen führt. Es zeigt 
sich leicht), daß dies nur für negative Werte von D der Fall; bei 
positivem D ergibt die Differentialgleichung bekanntlich immer, daß 
anfänglich schwankende Konzentrationen im Laufe der Zeit sich 
ausgleichen. 

Nun ergibt eine einfache thermodynamische Überlegung, daß 
für den Diffusionskoeffizienten D die Beziehung gilt: 


D=Cc 5° 


Nach der Diffusionstheorie ist nämlich die Zahl der Moleküle, die 
sich in der Zeiteinheit durch ein Konzentrationsgefälle hindurch 
über eine Strecke 42 hinweg bewegen, gleich: 


Andererseits ist entsprechend dem zweiten Hauptsatz diese Zahl: 
proportional der bei dem Vorgang eintretenden Abnahme des thermo- 
dynamischen Potentials; außerdem der Konzentration?), also gilt 


0 0G de 
x + Az 


de On Ox 


’ 


woraus die oben angegebene Beziehung für D folgt. 


Damit D oberhalb der Konzentration c = 4°/, negativ wird, 


muß also die Größe = bei dieser Konzentration ein Maximum haben 


1) Vgl. U. Dehlinger, a. a. O. Das von c abhängige D wird dort als 
„effektiver Diffusionskoeffizient“ bezeichnet. 


2) Dies gilt nur für Diffusionen, die auf "/orgängen 1.Ordnung beruhen. 
In Substitutionsmischkristallen ist statt dessen Proportionalität mit c(l—c) anzu- 
nehmen. Übrigens werden die Endergebnisse davon nicht berührt. Vgl. U. Deb- 
linger, Chemische Physik der Metalle und Legierungen, 1939. 
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und für ¢ > 4°/, mit wachsendem c abfallen. Dies läßt sich aber 
leicht verstehen: Es ist ja allgemein 


aG _ au 
On On 


eV, 
On 


6s 


Dabei ist — su gleich der Arbeit, die beim Aufbau der Lösung aus 


ihren einzelnen lonen und Molekülen nach außen abgegeben wird. 
Diese Größe wird in dem Konzentrationsgebiet nahezu von c unab- 
hängig sein, in dem jedes Ion noch von einer geschlossenen Wasser- 
hülle umgeben ist. Von c= 4°/, ab ist das nach dem obigen nicht 
mehr der Fall: mit wachsendem c werden jetzt in steigendem Maß 
je zwei Ionen unmittelbar benachbart sein, also starke Coulombsche 
Kräfte aufeinander ausüben, während vorher nur Dipolkräfte wirksam 


Abb. 4. Theoretischer Verlauf der 
Kurven, aus welchen sich das thermo- 
dynamische Potential der Lösung 
durch Addition zusammensetzt. Die 
Knickstelle entspricht c = 4°%,. 


— ist von T unabhängig 


waren. Von c = 4°/, ab wird also die Aufbauarbeit stark zunehmen, 


somit 5, abnehmen. Die Konzentrationsabhängigkeit der Größe 


-T a ist mit guter Näherung durch RT In c gegeben; p = kann 
vernachlässigt werden. In Abb. 4 ist se und RT In c aufgetragen. 
Man sieht, daß bei tieferen Temperaturen die Summe beider Aus- 
drücke ein Maximum für c = 4°/, besitzt, während bei hohen Tem- 
peraturen dieses Maximum allmählich nach c = 1 rücken wird. 
Wegen der statistischen Regellosigkeit der Verteilung von 
Wassermolekülen und Ionen in der Lösung werden auch schon unter- 


halb c= 4°/, einige Ionen sich nahe kommen. Das hat zur Folge, 
daß die Kurve für in Abb.4 keinen scharfen Knick bei c = 4°/, 


hat, sondern allmählich von einem Ast zum andern geht. Das Maximum 
der Summenkurve rückt dadurch zu etwas kleineren Konzentrationen 
und nähert sich dem Punkt c= 4°/, mit steigender Temperatur. 
Wir erhalten also für den Temperaturverlauf und die Lage des 
Maximums von . theoretisch gerade die Werte, die wir experi- 
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mentell für die Grenze der spontanen Keimbildung festgestellt haben, 
Damit dürfte der Zusammenhang dieser Grenze mit dem Negativ. 
werden des Diffusionskoeftizienten gesichert sein. 


Es ist noch zu bemerken, daß sich die Abnahme von = ober- 


halb c = 4°/, auch in der Löslichkeitskurve (Abb. 3) bemerkbar macht, 
Diese biegt in der Nähe dieser Kurve in dem zu erwartendem Sinn 
von der Kurvenform ab, die bei idealer Lösung, also von der Kon- 


zentration unabhängigem Zu erwarten wäre. 


Anschaulich ist also der Keimbildungsvorgang in der wäßrigen 
Lösung in folgender Weise zu beschreiben: Zunächst bilden sich 
an einzelnen Stellen unter dem Einfluß statistischer Schwankungen 
kleine Konzentrationshäufungen. ‘Unterhalb der Grenzkonzentration 
werden diese Anhäufungen durch die Temperaturbewegung immer 
wieder abgebaut. An der Grenzkonzentration selbst bleiben sie lange 
Zeit bestehen, da hier der Diffusionskoeffizient Null ist’). Oberhalb 
der Grenzkonzentration aber sind die Anhäufungen der Konzentration 
thermodynamisch stabil. Solange sie noch klein sind, ist ihr Sta- 
bilitätsgrad gering, so daß sie durch größere Schwankungen in der 
Geschwindigkeitsverteilung wieder zerstört werden können; einzelne 
bleiben aber übrig, wachsen im Laufe der Zeit weiter und werden 
dadurch stabiler. Da die Grenzkonzentration dadurch bestimmt ist, 
daß ein Teil der Ionen nur noch durch eine Schicht von Wasser- 
molekülen voneinander getrennt ist, können wir dieses Wachsen in 
folgender Weise atomistisch beschreiben. Aus den durch eine 
Wasserschicht getrennten Ionen entstehen durch Einzelsprünge Ionen- 
paare, die durch starke Coulombsche Kräfte zusammengehalten 
werden und daher verhältnismäßig lange Zeit bestehen können. Haben 
wir an einer Stelle eine Anhäufung von solchen Paaren, so werden 
aus den in ihrem Zentrum liegenden Paaren mit großer Wahrschein- 
lichkeit Gruppen von drei und mehr Ionen entstehen; die dadurch 
in der Nachbarschaft des Zentrums entfernten einzelnen Ionen werden 
durch Diffusion vom Rand der Anhäufung herangeschafft. 


VI. Schlußbemerkung 
Die vorliegende Untersuchung wurde ursprünglich durch folgende 
Fragestellung veranlaßt: Die fundamentalen chemischen Reaktionen 
der Biologie bestehen sicher darin, daß hochpolymere, kristallähnliche 


1) Die bekannten Dichtenschwankungen am kritischen Punkt entstehen 
thermodynamisch ganz ähnlich wie diese stabilen Konzentrationsschwankungen. 
Auch bei unserer Grenzkonzentration hört die „thermodynamische Stabilität“ 
der Lösung auf. 
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Stoffe, z. B. die Gene im Zellkern oder die Viren, sich katalytisch 
oder autokatalytisch aus kleineren Bestandteilen aufbauen, die in 
wäßriger Lösung vorliegen. In den letzten Jahren wurde nun mehr- 
fach die Frage diskutiert, ob die dabei wirksamen Kräfte ganz spe- 
zifische seien, z.B. außergewöhnliche quantenmechanische Resonanz- 
kräfte zwischen den Molekülen, oder ob es sich um normale Atom- 
kräfte z.B. Dipolwirkung handle, wie sie auch beim Aufwachsen von 
sonstigen Kristallen wirksam sind!) Gegen die letztere Auffassung 
war einzuwenden, daß sie nicht begründen könne, warum bei den 
biologischen Reaktionen niemals eine spontane, nichtkatalytische 
Bildung der genannten Stoffe eintritt. Die vorliegende Untersuchung 
zeigt nun, daß die Konzentration der spontanen Keimbildungsgrenze 
dadurch bestimmt ist, daß jedes der in Frage kommenden Komplexe 
gerade noch eine vollständige Wasserhülle besitzt. Bei großen Kom- 
plexen, um die es sich bei organischen Stoffen stets handelt, liegt 
diese Konzentration außerordenlich hoch, so daß sie unter bio- 
logischen Verhältnissen überhaupt nicht mehr vorkommt. Wir 
haben also bei Lösungen größerer organischer Moleküle oder Radi- 
kale überhaupt nicht mit spontaner Keimbildung zu rechnen); der 
oben erwähnte Einwand gegen die Annahme normaler chemischer 
Kräfte bei biologischen Aufbauprozessen entfällt somit. 


Zusammenfassung 


1. In einer durch mehrmaliges Umkristallisieren systematisch 
gereinigten Alaunlösung werden bei 0°, 20° und 40° C die Warte- 
zeiten bestimmt, nach welchen spontane Keimbildung beginnt. Es 
zeigt sich, daß unterhalb einer wenig von der Temperatur ab- 
hängenden, nahe bei 4 Mol-Prozent KAI(SO,), gelegenen Konzentrations- 
grenze auch nach mehr als einem Jahr keine spontane Keimbildung 
eintritt. Es ergibt sich also eine scharfe Grenze zwischen meta- 
stabilen und stabilen übersättigten Lösungen. 

2. Ein Versuch, die für Keimbildung aus dem Dampfzustand 
geltende Theorie von Volmer u. a. auf diese Grenzkurve anzuwenden, 
zeigt, daB sowohl die Schärfe der experimentellen Grenze als auch 
ihr Temperaturverlauf mit dieser Theorie nicht zu vereinbaren ist. 


1) Vgl. U. Dehlinger, Naturwiss. 23. S. 558. 1935; 25. S. 138. 1937. 
K.Sommermeyer und U. Dehlinger, Phys. Ztschr. 40. S. 67. 1939. 
Th. Neugebauer, Phys. Ztschr. 40. S. 406. 1939. 

2) In einzelnen Fällen ist dies auch experimentell bestätigt. So zeigt nach 
F.8.Kipping und W.J. Pope, J. chem. Soe. 95. 103. 1909 Natrium-Ammonium- 
Tartrat auch nach Einengung bis zu Sirupdicke keine spontane Keimbildung. 
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3. Dagegen kann die zuerst bei der Ausscheidung in metallischen 
Mischkristallen aufgefundene Keimbildung durch negative Diffusion 
dies ohne Schwierigkeit erklären. Man muß dazu annehmen, daß 
jedes gelöste Ion sich mit 6 Wassermolekülen zu umgeben sucht; 
oberhalb 4 Mol-Prozent ist dann ein Teil der Ionen nur noch durch eine 
Schicht von Wassermolekiilen getrennt. Die dann auftretenden 
Coulombschen Kräfte zwischen den Ionen verursachen bei tieferen 
Temperaturen ein bei 4 Mol-Prozent liegendes Maximum im Konzen- 
trationsverlauf des thermodynamischen Potentials und damit einen 
negativen Diffussionskoeffizienten oberhalb dieser Konzentration. 


Stuttgart. 
(Eingegangen 29. Januar 1941) 
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